CERN-PH-TH/2007-260 



Magnetized CMB observables: a dedicated numerical approach 



Massimo Giovannini'^''^ and Kerstin E. Kunze'*''^ 



""Centra "Enrico Fermi" , Compendia del Viminale, Via Panisperna 89/A, 00184 Rome, Italy 

* Departamento de Fisica Fundamental, 
Universidad de Salamanca, Plaza de la Merced s/n, E-37008 Salamanca, Spain 

Department of Physics, Theory Division, CERN, 1211 Geneva 23, Switzerland 



Abstract 

Large-scale magnetic fields affect the scalar modes of the geometry whose ultimate effect is 
to determine the anisotropies of the Cosmic Microwave Background (CMB in what follows). For 
the first time, a consistent numerical approach to the magnetized CMB anisotropies is pursued 
with the aim of assessing the angular power spectra of temperature and polarization when the 
scalar modes of the geometry and a stochastic background of inhomogeneous magnetic fields are 
simultaneously present in the plasma. The effects related to the magnetized nature of the plasma 
are taken into account both at the level of the dynamical equations and at the level of the initial 
conditions of the Einstein-Boltzmann hierarchy. The temperature and polarization observables are 
exploited to infer the peculiar signatures of a pre-equality magnetic field. Using the extrapolated 
best fit to the three year WMAP data the increase and distortions of the first seven peaks in the 
TT autocorrelations are monitored for different values of the regularized magnetic field intensity 
and for the physical range of spectral indices. Similar analyses are also conducted for the first few 
anticorrelation (and corrrelation) peaks of the TE power spectra. Possible interesting degeneracies 
and stimulating perspectives are pointed out and explored. 



1 The general framework 



The origin of large-scale magnetic fields is still a poorly understood subject as much as the origin of 
the primordial spectrum of the density contrast was unclear almost two score years ago when the first 
attempts of understanding the origin of large-scale structure actually begun. The first observational 
evidence of the existence of large-scale magnetic fields in our galaxy goes back to the work of Hiltner 
[1] and Hall [2] (see also [3]). The problem, in itself, has many facets which can be summarized as 
follows: 

• observational evidence of the existence of large-scale magnetic fields is indisputable from a number 
of different observations ranging from Faraday rotation to synchrotron emission; 

• objects of different sizes (and different evolutionary histories) possess sizable magnetic fields in 
the /i G range: a striking example of this statement are, for instance, spiral galaxies, elliptical 
galaxies and a class of x-ray bright Abell clusters; 

• theoretical evidence of the mechanism responsible of the degree of magnetization of the present 
Universe is still under active discussion. 

We have a rather plausible control of the dynamics of electromagnetic fields in plasmas: since the pi- 
oneering work of Alfven [4] the various descriptions of weakly coupled plasmas have been tested both 
in astrophysical systems and, most importantly, laboratory experiments [HI El [TJ [8] . Various plasma 
descriptions (covering different branches of the spectrum of plasma excitations) allow to predict insta- 
bilities in terrestrial tokamaks. An example of successful (but not unique) framework is the one-fluid 
plasma theory which, under some circumstances, reduces to the well known magnetohrdrodynamics 
(MHD in what follows). All the astrophysical attempts for the justification of large scale magnetic 
fields rest on the assumption that, at some time prior to galaxy formation, appropriate initial con- 
ditions for the MHD evolution should be present. While this point of view is shared by the whole 
community, opinions vary as far as the real primordial nature of magnetic fields is concerned. For 
informative reviews on this broad subject the interested reader may usefully consult the publications 
reported in [9] (see also [TOj for more observational aspects of large-scale magnetization). 

Theoretical diatribes cannot decide wether or not large-scale magnetic fields are a cosmic relic. It is 
more rewarding, in our opinion, to apply Occam's razor and pose a more modest but yet experimentally 
answerable question: were large-scale magnetic fields present prior to matter-radiation equality? The 
latter question is answerable since, as it will be specifically shown, the effects of large-scale magnetic 
fields can be read-off from the properties of CMB observables. 

In recent years diverse data sets seem to conspire towards a sort of paradigm which is customarily 
employed in the interpretation of cosmological data. This paradigm is often dubbed as ACDM lore 
where A stands for the dark-energy component (parametrized as a cosmological constant) and CDM 
stands for the cold dark matter component. The data sets supporting this general view range from 
the CMB data0 to the large-scale structure data [TTl [TS] and, finally, to the type la supernova data in 
their various incarnations p^l20] (see also [21] for weak lensing data). The combined analysis of these 
three classes of data support the ACDM paradigm. A primordial spectrum of adiabatic perturbations, 

^ Among them, the WMAP data [TTJ [T^l [131 HH US] , the data of other balloon borne experiments as well as the data 
of various terrestrial arrays [16] . 
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present for typical wavelengths larger than the Hubble radius prior to Hydrogen recombination and 
after matter-radiation equality, is the main responsible for the normalization of the CMB temperature 
autocorrelations (TT correlations in what follows). The cross-correlations between temperature and 
polarization (TE correlations in what follows) lead to a typical anticorrelation peak [Hj which is 
often quoted as the golden signature of the predominant adiabaticity of the CMB initial conditions. 
The polarization autocorrelations (EE correlations in what follows) have been partially observed by 
the WMAP experiment itself [12] as well as by other terrestrial arrays. The viability of the ACDM 
paradigm is customarily completed by the addition of various parameters whose explicit determination 
can either confirm or improve the pure ACDM scenario. Among these models we can mention, just as 
an example, the addition of a tensor component, the inclusion of a specific form of the barotropic index 
for the dark energy component (different from the one of a pure cosmological constant), the search 
for specific corrections in the scalar power spectrum motivated, with various degrees of theoretical 
accuracy, by the (yet unknown) physics of the Planck energy scale (see, for instance, some of the 
analyses reported in fl]). There is the hope that the unprecedented accuracy of the data of the 
Planck explorer mission [22] will allow to rule in (or rule out) some of these theoretical possibilities B 

The aim of this paper is to complement the ACDM paradigm with the presence of a pre-decoupling 
magnetic field (see also [25]). Such a completion is not only motivated theoretically but also observa- 
tionally: since we do observe large-scale magnetic fields later on, it is plausible to posit, as falsifiable 
hypothesis, their existence also prior to equality. Definite answers to such a question will come, in 
this context, by confronting the completion of the ACDM paradigm with all the available cosmological 
data in the same way as, for instance, a pristine tensor contribution to the CMB anisotropies can be 
constrained by adding, in the parameter estimation, a stochastic background of long-wavelength gravi- 
tons. In this investigation one of the first goals of the program will be reached: a dedicated numerical 
approach for the magnetized CMB observables will be constructed. Without this step sound strategies 
of parameter estimation will be forlorn. This step is often rather straightforward in different cases 
since all the available codes contemplate, for instance, the inclusion of tensor modes or the inclusion 
of a peculiar barotropic index for the dark energy component. However the latter statement does not 
apply to the case of large-scale magnetic fields. 

The analysis of the interplay between large-scale magnetic fields and CMB observables might be 
traced back to the pioneering works of Zeldovich [26] and Harrison [27]. In recent years it has been 
understood that large-scale magnetic fields may affect the vector and tensor modes of the geometry 
and may also affect, indirectly, the CMB polarization [2S1 [2H] (see, for instance, [30] for a topical 
review on this subject). The main obstacle to a systematic analysis of the current data in the light of 
a magnetized component is represented by our lack of understanding of the close relationship between 
the large-scale magnetic fields and the scalar modes of the geometry. Indeed, as explicitly suggested 
by observations, the inhomogeneities in the CMB temperature and polarization can be attributed to 
curvature perturbations. 

An impasse then seems to arise. The vector and tensor modes induced by large-scale magnetic 
fields are very small at large length scales (i.e. small multipoles). It is thus rather hard to imagine the 
possibility of including the magnetic field contribution as an appropriate fit parameter in an extended 
version of the ACDM paradigm if the effects of magnetic fields on the scalar modes of the geometry 
are unknown. 

^Various models going beyond the pure ACDM scenario can be confronted with the foreseeable sensitivities of the 
high and low frequency instruments embarked on the Planck spacecraft. Sec, for instance, [23[ I24j . 
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A theoretical framework for the analysis of scalar modes of the geometry and large-scale magnetic 
fields has been recently developed in a series of papers ranging from the accurate determinations of 
adiabatic initial conditions in the presence of a scalar mode of the geometry [31] , to the extension of the 
tight- coupling expansion [32] and to the semi-analytical calculation of temperature autocorrelations 
[33] . One of the aims of the present paper is to translate into a dedicated numerical approach all the 
theoretical understanding of the interplay between the scalar modes of the geometry and the large scale 
magnetic fields that has been pursued, through various steps, in [311 \32\ 133] . The CMB effects related 
to the scalar modes of the geometry are the most difficult ones already in the absence of large-scale 
magnetic fields. The complication comes from the physical observation that the density contrasts of 
the various species do not couple to the tensor and the vector modes but they do couple to the scalar 
modes of the geometry. Conversely the curvature perturbations are the source of the evolution of 
the density contrasts for photons, for baryons and for the CDM species. In the magnetized case this 
occurrence is even more acute since large-scale magnetic fields and plasma effects propagate both at 
the level of the Boltzmann hierarchy and at the level of the perturbed Einstein equations. 

The approach studied in the present investigation is based on a faithful MHD description of the pre- 
decoupling plasma. More refined descriptions of the high frequency branch of the spectrum of plasma 
excitations (valid for frequencies comparable with the electron plasma frequency) must reproduce any- 
way the approach described here when the typical length-scales and dynamical times are, respectively, 
much larger than the Debye length and much larger than the inverse of the plasma frequency [HI [6] . 

The plasma will not only be populated by an electromagnetic component but also by fiuctuations 
of the geometry which should be treated relativistically since their typical wavelengths, at the onset 
of the numerical calculation, will be much larger than the Hubble radius at the corresponding epoch. 
The description of the scalar modes of the geometry will be formulated in the peculiar language of the 
synchronous coordinate system. One of the first codes developed for the analysis of CMB anisotropics is 
COSMICS [311 [35]. The CMBFAST [361 [37] code is based originally on COSMICS and in many respects 
it can be said that COSMICS is the ancestor of CMBFAST. As in COSMICS, also in CMBFAST the 
dynamical evolution across equality and decoupling is followed in the synchronous coordinate system. 
The synchronous approach carries necessarily a possible ambiguity on the complete removal of the 
gauge freedom. It is actually known since the early eighties [38] that, in the synchronous gauge, the 
coordinate system is only fixed up to a pair of space-dependent integration constants. The remaining 
gauge freedom must be removed from the initial conditions of the Boltzmann hierarchy to avoid the 
dangerous presence of spurious (i.e. gauge) modes. The way to handle this potentially annoying 
problem resides in the abihty of treating the same problem also in contiguous gauges where the gauge 
parameters are completely fixed. An example of this technique is the longitudinal gauge [271 ESI ES] 
(see also [40]) which has been also exploited for this purpose in the absence of magnetic fields. This will 
also be the path followed in the present analysis (see also [SI1132]) by including consistently the effects 
due to the global magnetization of the plasma. Finally, in the present study the scalar vector and 
tensor modes are separated on the basis of their transformation properties under three-dimensional 
rotations jEH] (see also [H]). There exist approaches which are fully covariant [12] and which have been 
also applied to the case of large-scale magnetic fields [13] without leading, however, to any explicit 
estimate either of the Sachs- Wolfe plateau or of the temperature autocorrelations in the Doppler region 
as reported in pTIJ. 

With these specifications in mind, our code is based on CMBFAST with integration along the 
line of sight. The main difference stems from the consistent introduction of a magnetized component 
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both al the level of the initial conditions and at the level of the dynamical equations. This choice has 
been also dictated by the fact that the WMAP collaboration used also CMBFAST for the analysis 
of the observational data. It is plausible to think that this numerical approach can be also effective 
throughout the next decade for the analysis of forthcoming data like the ones of the Planck explorer 
[22]. 

The plan of our paper is therefore the following. In Section [2] we will review the main evolution 
equations to be integrated. Particular attention will be given to the way large-scale magnetic fields are 
included in the pre-equality and pre-decoupling physics. In Section [3] a theory of the magnetized initial 
conditions will be formulated when the dominant source of curvature inhomogeneities is the standard 
adiabatic mode. Section H] is devoted to the calculation of temperature autocorrelations. The results 
on the polarization observables are collected in Section El The distinctive features introduced in the 
angular power spectra by the presence of large-scale magnetic fields will be further scrutinized in 
Section [6l Section [7] contains our concluding remarks. In the Appendix we collected auxiliary material 
on the longitudinal gauge description sticking, however, only to those themes that are germane to our 
calculations. 



2 The full content of the magnetized plasma 

Below the temperature of neutrino decoupling, the content of the plasma is formed both by neutral 
species and charged species. The neutral species are cold dark matter (CDM in what follows), neutrinos 
(which will be taken to be massless) and photons. The charged species are baryons and electrons. The 
three observational data sets (i.e. CMB anisotropies [HI [T2], large-scale structure [T71 HB] and type 
la supernovae [IHIEQ]) suggest, that a cosmological constant term should also be added. The pivotal 
ACDM lore is then delicately improved by positing, already prior to equality, that the primeval plasma 
is effectively magnetized. 

Neutral species are indirectly affected by the presence of large-scale magnetic fields. In fact, 
magnetic fields gravitate and contribute both to the Hamiltonian and momentum constraints as well as 
to the dynamical evolution of the gravitational inhomogeneities. The evolution of the charged species 
will be followed using MHD [51 El E] generalized to the situation where the geometry is dynamical and 
where gravitational inhomogeneities are also simultaneously present. The adoption of this scheme is 
dictated by the typical hierarchy of the Coulomb and Thompson scatterings. 

The position of the first Doppler peak implies, in a ACDM framework, that the background geom- 
etry is spatially fiat. The inclusion of spatial curvature amounts to an overall shift of the TT power 
spectra and hence to a change in the position of the first peak. The line element can then be written, 
in the conformal time coordinate r, as: 

a' 

ds^ = a^{T)\dT^ - dx\ a{T)dT = dt, H = - = aH. (2.1) 

a 

In Eq. (12. ip . t denotes the cosmic time coordinate and H the Hubble rate; the prime will denote, 
throughout the paper, a derivation with respect to the conformal time coordinate r. The evolution of 
the scale factor a(r) is determined by the appropriate Friedmann-Lemaitre equations: 

= -—a^pt, (2.2) 
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-n' = ATiGa\p, + p,), (2.3) 
p; + 37^(pt+Pt) = 0, (2.4) 

where pt and pt denote, respectively, the total energy density and pressure of the plasma, i.e. 

Pt = pc + Pu + + Pe + ph + Pk, Pt = y + y-pA- (2.5) 

The subscripts in Eq. (12.51) refer to the various components of the plasma mentioned in the first 
paragraph of the current section. 

The time scales involved in the present study are the ones encountered in CMB physics: the equal- 
ity time (at which the radiation and the matter component have equal weight in Eq. (12. 2p ). the 
recombination time (at which the ionization fraction drops because neutral Hydrogen is formed), the 
decoupling time (at which the mean free path of the photons becomes comparable with the Hubble 
radius). The exact solution of Eqs. (12.21) . (12. Sp and (12. 4p in the absence of dark energy component 
(which can be neglected around decoupling) stipulates that the scale factor interpolates, in the confor- 
mal time coordinate r, between a linear evolution (typical of the pre-equality plasma) and a quadratic 
evolution (typical of the plasma around recombination and decoupling): 



^1 = 7^^7^- 288.15 (^)"'mpc. (2.6) 




In Eq. (12.61) cjm = ^o^mo where /iq (of the order of 0.7 in the standard ACDM framework) is the 
current indetermination on the Hubble rate. Given a generic species X, ujx = h'^^xo'- while VLxo is 
proportional to /ig ^) is, by its definition, independent of ho (this is the reason why it is sometimes 
called Hubble-free critical fraction). To estimate ti in Eq. (12. 6p it has been used that u;^ = 2.47 x 10^^, 
= 1.68 X 10-5 and = + u.^ = 4.15 x 10"^ In the code Eqs. (E3), O and ([23D are 
integrated numerically. Equation (12.61) can be anyway used for semi-analyitical estimates [22] and for 
the important problem of correctly setting the initial conditions of the Einstein-Boltzmann hierarchy 
(see Section[3]). According to Eq. (12. 6p Teq = (y^— l)ri ~ ti/2. The redshift is defined as 1 + z = a^/a 
by fixing = 1- The recombination redshift (taken for instance between 1050 and 1100) will determine, 

via Eq. (12. 6P Tree (note that Tree > Teq)- 

The coupled evolution of the fluctuations of the geometry (12. ip with the fluctuations of the plasma 
quantities will determine, thanks to the interaction with the magnetic fields, the peculiar features 
of the magnetized CMB anisotropies. The synchronous coordinate system is most easily presented 
already in Fourier space and it can be written a^ 



4^,(/c, r) = a\r) hk,h[K r) + G^k, r) [ kik, - l^., 



(2.7) 



where ki = ki/\k\ denotes the direction of the Fourier wave- vector and 6s reminds that we are consid- 
ering here only the effects of the scalar modes of the geometry which are, as already mentioned, the 
most difficult (but also the most relevant) when the plasma is effectively magnetized. 



■^Our conventions on the perturbations are summarized by Eq. (|2.7p . It should be borne in mind that the signature 
of the metric is mostly minus (see Eq. (12.11) ). Various treatments of this problem adopt the opposite convention (i.e. 
the signature (— , +, +, +)). 



5 



2.1 CDM particles and neutrinos 



Defining as 6c and 6^ the density contrast of CDM particles and tlie tliree-divergence of the corre- 
sponding pecuhar velocity, the evolution equations of the CDM sector are: 

S'c = -Oc + j, 9'^ + n9c = 0. (2.8) 

In spite of the fact that the CDM velocity field will be used, in Section [3] to remove partially the 
remaining gauge freedom, it is relevant to appreciate that Eq. (12. 8p can be also written as 

5c-|)" + h(5c-|)' = 0. (2.9) 

Defining, in analog terms, 6i, and 6i, as the neutrino density contrast and as the three-divergence of 
the neutrino peculiar velocity, the corresponding evolution equations are: 

K = -lOu + lh', (2.10) 
C = -A;V, + y(5„ (2.11) 

where a^, is the neutrino anisotropic stress (also appearing in the perturbed Einstein equations) which 
is related to the quadrupole of the (perturbed) phase space distribution as a^, = Ty^l'^'^ is the 
octupole of the (perturbed) phase space distribution. The presence of the quadrupole and octupole 
reflects the occurrence that neutrinos are coUisionless below temperatures of the order of the MeV 
and should therefore be treated in the framework of the appropriate Boltzmann hierarchy. Equations 
fl2.10p . (12. lip and (12.121) couple together the lowest multipoles and will be the ones used, in Section 
[3l to set initial conditions of the CMB anisotropics in the pre-equality regime. At later time, in the 
code, the neutrinos will be integrated using the perturbed form of the coUisionless Boltzmann equation 
written in the synchronous coordinate system: 

,2 



-r + y(/^' + 60 



(2.13) 



where fi = k ■ n. Equations (I2.10p . (12.111) and (I2.12p can be derived from Eq. (12.131) by expanding JF^ 
in series of Legendre polynomials according to the convention: 

^,(fc, n, r) = E(-^)'(2^ + ^)J'Ak, r)Pi{fi). (2.14) 

e 

For larger multipoles Eqs. (I2.10p . (12.111) and (12.121) can be written as 

Ki = ^^[^-^K^-D - + ^ > 3. (2.15) 

Equation (12.151) can be derived from Eq. (I2.13P if we multiply both sides by -P^(yu), i.e. by a generic 
Legendre Polynomial. By then integrating over fi (between —1 and 1) and by using the convention 
established in Eq. (12.140 . Eq. (I2.15P follows by appropriate use of the recurrence relation of the 
Legendre polynomials stipulating that {£ + l)P^+i(/i) = {2i + l)fiPe{fi) - iPi-iifJ.) US 09]. The 
numerical integration will demand to cut the hierarchy at an appropriately (large) multipole. The 
neutrino fraction in the pre-equality plasma will be denoted by and it is defined as: 

^. = T— , r = -Nj-) =0.681 (2.16) 



1 + r' 8 viiy 

By definition, the photon fraction at the corresponding time will be given by i?^ = 1 ~ Ru. 
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2.2 Pre-decoupling plasma 



Overall, the plasma will obey the Gauss constraint whose explicit form can be expressed, in the present 
case, as 

V ■ E = Aneirii -He), E = a^S, (2.17) 

where rig is the electron concentration and is the concentration of the ions (to be soon identified with 
protons); E is the electric field rescaled through the second power of the scale factor 0. Electrons and 
ions are approximately in thermal equilibrium for temperatures smaller than the MeV. The electric 
field appearing in the Gauss constraint will be screened for length scales L > Ad: 



AD(re,) = J-^ = 2.8f^rY^r''V^r m, (2.18) 
""^^ V 87re%o VO.I34; V0.023y VeVy ' ^ ' 

where T ~ Tc ~ Ti ^ T^. In Eq. f l2.18p the common value of the electron and ion concentrations is 
simply given by 

no = ne = ni = r]^n^, r/b = 6.27 x 10"^° {^^^ > (2-19) 



.0.023. 

where the ionization fraction is set to 1 as it is consistent for pre-recombination temperatures. Since 



the Hubble rate at equality is given by i^cq = 1-65 x 10 ^^(u;m/0.134)2Mp we have that 



HMr^a) ^ 2.85 X lO-^^ ( ( {—\\ (2.20) 
''''^ V0.134/ V0.023/ VeV/ ^ ' 

The hierarchy between the Debye length Ad and the other length-scales of the problem persists if we 
move from the equality time to the decoupling time. According to Eq. (I2.18p Ad scales § as a(r) (since 
T ~ and — a~^^. Thus A(rdcc) — 3(co'M/0.134)Aeq where the factor 3(a;M/0.134) arises because 
(odec/oeq) — 2.393(071^1/0.134) (taking, as an example, l-|-2;dec — HOO and fixing the ionization fraction 
as Xe — 1). 

Recombination entails a sudden drop in the ionization fraction. From the usual considerations 
involving Saha's equation, around decoupling, ~ 10"^. The Debye length increases then (see Eq. 
f l2.18p ) by a factor 10^'^ which is still minute in comparison with all the other lengths of the problem. 
Note that AD(Tdcc) is not only parametrically smaller than the Hubble radius, but it is also negligible 
in comparison with the sound horizon at the corresponding epoch, i.e. 

rs(rdec) = / c?rc,b(r) = / , (2.21) 

/3[l + i?b(r)] 



where Csb(T) (the characteristic sound speed of the baryon-photon system in the tight coupling ap- 
proximation) is defined in terms of the baryon to photon ratio i?b (see also Eqs. fl2.83p and fl3.ll) ): 

i?b(.) = ^^= 0.664f^Vl^V (2.22) 
^ ^ 4p^ V0.023A2 + 1/ ^ ' 



^In the present paper we will denote with calligraphic style the electric and magnetic fields which are not rescaled 
through 0? , i.e. E = a^E and B = a^B. 

""In principle the electron and baryon temperature (in the absence of energy exchange and in the absence of photons) 
should scale as a~^. However, as it will be discussed in a moment, the Coulomb and Thompson rates of interactions are 
both larger than the Hubble rate and Tc ^ ~ — T (see Eqs. (|2.30l) - (|2.35|) and discussions therein). 
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Since the dark energy component is negligible around decoupling the integral appearing in Eq. (12.211) 
can be estimated analytically and the overall result can be expressed as: 

rs(rdec) 2998 2 



In 



a/1 + Ci + v^Ci + C1C2 



1 + yclc^ 



(2.23) 



Mpc VI + -2dec V3 t^MCl 

where 

07 K / 1100 \ 0.045 /l + 2d^\ 

With our fiducial values of the parameters, rs{T^ec) lies between 150 and 200 Mpc. But now this figure 
should be compared with the Debye length A(rdcc) — 2.5 x 10'^ m (having taken into account the drop 
in the ionization fraction). Thus, as anticipated, AD(Tdcc)/''"s(''dcc) — 10~^^. So the plasma is, to a very 
good approximation globally neutral. 

As already mentioned, for temperatures smaller than the temperature of neutrino decoupling 
baryons and electrons interact strongly through Coulomb scattering. The corresponding rate, for 
T > eV, is0: 

^i^^4xlO^^Xef-r''V^V (2.25) 
H \eVj V0.023y ^ ^ 

Since after equality H oc T^/^ ( and ignoring for the moment the drop in the ionization fraction) the 

ratio of Eq. (I2.25P gets frozen. Equation (I2.25P justifies to consider a unique baryon-lepton fiuid as a 

single dynamical entity. In Eq. (I2.25P the Coulomb rate has been computed by recalling that 

Tcoul = ''^ea^e'^thC^Coul, 

.« = -flnAe. Ae = 5(-) ^. (2.26) 

where In Ac — 14.71 for typical values of cUb- The Coulomb cross section is the main responsible for 
the conductivity of the plasma which can be estimated as the ratio between the square of the plasma 
frequency and the Coulomb rate (which is also, by definition, the collision frequency), namely 

^ciT) = -^ = —^T(-L.) ^0.16(4-) eV, (2.27) 
47rrcoui aemlnAc KnieJ VeV/ 

where, as it should, the electron concentration effectively simplifies in the final expression. It is useful 
also to estimate, at this point, the plasma frequency of the electrons, i.e. 




o^pe = = 28.05Vi;f-W^V'' MHz. (2.28) 

V "^e VeV/ V0.023/ ^ ' 

Thus the typical length-scales are much larger than the Debye scale. The typical time-scales greatly 
exceed uj~}. This is the realm of MHD. 

Prior to equality electrons and protons interact also with photons via Thompson cross-section. Pro- 
tons can be neglected in the Thompson mean free path (determined by electron-photon interactions); 
the Thompson rate in units of the Hubble rate is, for T > eV, 



H 



5.9x lO^Xef^^V^V (2-29) 
V0.023/VeVy ^ ^ 



^In the case of a proton (or of an electron) impinging on an electron (or on a proton) the Rutherford cross section is 
logarithmically divergent at large impact parameters when the particles are free. In the plasma around decoupling the 
logarithmic divergence is avoided because of the Debye screening length: the cross section is then known as Coulomb 
cross section and the logarithmic divergence is replaced by the so-called Coulomb logarithm. 
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After equality Tt^^/H is proportional to T^/^ and it becomes eventually much smaller than one as 
Xe drops at recombination. Deep in the radiation epoch, i.e. when the initial conditions of CMB 
anisotropies are set numerically, the Coulomb rate of Eq. (12.25^ is larger than the Thompson rate but 
while the Thompson rate increases with the temperature (see Eq. fl2.29p ) the Coulomb rate decreases. 
The meeting point of the two rates occurs close to the MeV. Initial conditions will then be set in the 
radiation epoch when both Coulomb and Thompson scattering are large. Also photons are strongly 
coupled to the baryon-lepton fluid. So, a unique physical entity emerges, i.e. the so-called baryon- 
lepton-photon fluid. This fluid is often dubbed as the baryon-photon fluid by implicitly including the 
electrons in the baryonic component thanks to the strength of Coulomb coupling. 

The various species have all a putative common temperature T, i.e. Te ^ Tp ~ = T. This 
statement will now be justifled. Electrons and protons, being massive, have energy densities and 
pressures which can be written, respectively, as 



Pe = rie 
Pp = np 



rrie + -Te 



Pe = rieTe, (2.30) 
Pp = UpTp. (2.31) 



The photon energy density and pressure will be instead, as it is well known, = (7r^/15)T^ and 
p^ = p^/3. Ignoring, for the moment, the other species of the plasma, the total conservation equation 
assumes the form d{a^p)+pd{a^) = where p = {pc+ Pp+ P'y) andp = (Pe+Pp+P^y)- Since = Up = uq 
(and both scale as a~^) the total conservation equation implies the following differential relation: 

aXd{aT^) + d[a\n + Tp)]=0, A = 2.8 x 10^ (^^^) (2.32) 

Up to numerical factors A is the ratio between the entropy density of the photons and uq which is, in 
turn, roughly 10 orders of magnitude smaller than the photon concentration. The electron and proton 
temperatures vary adiabatically as (in the absence of photons) while the radiation temperature 
varies as (in the absence of protons and electrons). But we do know that electron-photon and 
electron-proton interactions tie the temperatures close together. Equation (12.321) can then be solved 
assuming, to lowest order, Tg ~ Tp ~ = T. The differences in the various temperatures can be 
estimated: if the differences are small the assumption of a common temperature is justifled. To lowest 
order Eq. (12.321) becomes 

cilnT^ A + 4 ^ -, 2 , , 

——^ = -- , T,^a-^-A 2.33 

d\na A + 2' f ^ v ; 

where the second relation is obtained from the flrst one after expanding the obtained result in powers 

of 1/A. Equation (12.331) shows that, indeed, Ty evolves in a way which is intermediate between 

(as implied in the absence of electrons and protons) and a~^. The differences between the T^, Tg and 

Tp depend upon the Thompson and Coulomb rates. The rate of gain of energy per electron as well as 

the rates of gain of electron and proton thermal energies can be written as 

1 d{a^Te) 



T, — ~rTh(Te — Tj) — FcouK^e " ^p), 

dt 

ld{a%)_ ^ .r^ r^, Xd{aT^^ 



a"^ dt a dt 

Equation (12.341) implies, as expected, that 



rcoui(Tp - Tg), '—^ = -rTh(T, - Tg), . (2.34) 



T^ — Tp H Tp — Tr, H 

^—^^2 — , (2.35) 

J- J- Th J- J- Coul 
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The same hierarchy between Coulomb and Thompson scattering rates also determines the small tem- 
perature differences between electrons, protons and photons. Thus electron-proton collisions are suf- 
ficiently fast to assess that electrons and protons have indeed the same putative temperature. 



2.3 The baryon-lepton-photon fluid 

The strength of the Coulomb coupling implies that, effectively, there is a unique velocity field which 
is the centre of mass velocity of the electron-proton fluid, i.e. 

Vh = : • (2.36) 

rrte + rup 

The velocity is the bulk velocity of the plasma O [6]. The evolution equation of vy, can be obtained 
by summing up the evolution equations of electrons and ions as it happens in the usual MHD 
treatment of the problem. In the synchronous gauge, the baryon velocity and the baryon density 
contrast obey, respectively, the following pair of equations 



5'. = -0, + ^ + (2.37) 
2 a^ph 



e', + ne, = ^^e'ie, - e,) + Ydl^, (2.38) 

3 Pb a^Pb 

where, the divergence of the baryon velocity and the differential optical depth have been introduced 
as 

= V ■ Vh, e' = Xe (TT'^e- (2.39) 

ao 

Equation fl2.38l) has been written, unlike the analog equations for CDM and neutrinos, not in Fourier 
space but in real space to emphasize the presence of a new term which is nothing but the MHD 
form of the Lorentz force given hj J x B where J is the Ohmic current and B is the magnetic field. 
The electric field appearing in Eq. (]2.37p is negligible since, in the plasma frame, the conductivity 
effectively suppresses the Ohmic electric fields. 

According to Eq. fl2.38p the baryon-lepton fluid exchanges momentum with the photons. The 
lowest two multipoles of the Boltzmann hierarchy of the photons, namely the density contrast (i.e. the 
monopole) and the three-divergence of the velocity field (related to the dipole of the intensity of the 
brightness perturbations) are: 

6' = -\e^ + lh\ (2.40) 



3 ' 3 
1, 
4 



^; = --wH, + e\eh-e,). (2.41) 



While the sum of the electron and proton equations leads to Eq. (12.381) . their difference leads to the 
Ohm law [5^, 6] which relates the total current J to the electric field through the conductivity, i.e. 



J = a{E + VhX B), a = a^a, J = a-^j, (2.42) 

where (Tc denotes the fiat-space conductivity; furthermore, as already mentioned, B = a^B and E = 
a?£. The usefulness of the latter rescalings can be understood by looking also at the other MHD 
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equations, namely 



V ■ ^ = 0, V ■ = 0, (2.43) 

V X B = Air J, — + V X E = 0. (2.44) 

or 

Equations fl2.43p and (12.441) have the same form they would have in flat space. The space-time is 
however curved and with line element given by Eq. (12. ip . Maxwell equations in conformally flat 
backgrounds are known to be invariant under a Weyl rescaling of the metric. Consequently the 
corresponding evolution equations have exactly the same form they would have in Minkowskian space- 
time provided the fleld are appropriately rescaled and provided the conformal time coordinate r is 
consistently employed. 

The displacement current does not appear in Eq. (I2.44p . Indeed MHD is a description that holds 
for typical length-scales that are larger than the Debye length and for typical time-scales that are 
much larger than the inverse of the plasma frequency. In other words, if we are interested to study 
the high frequency branch of the spectrum of plasma excitations we should resort to a full kinetic 
(Vlasov-Landau) description [B]. The Ohmic current can then be related to the magnetic fleld, i.e. 

J=-^Vx5, V-J = 0. (2.45) 

An 

Thus, the total current, the electric fleld and the magnetic fleld are all solenoidal albeit for rather 
different physical reasons. Equation (I2.45P can be used to compute explicitly the Ohmic electric fleld, 
i.e. 

E = -v^x B + ^-^, (2.46) 

4770" 

which shows that electric flelds vanish, at flnite conductivity, in the baryon rest frame. They are 
therefore smaller than the magnetic flelds since, as previously shown explicitly, the pre-decoupling 
plasma is an excellent conductor. The latter statement deflnes the plasma frame, i.e. the frame where, 
thanks to the large value of the conductivity, the electric flelds vanish while the magnetic fleld are not 
dissipated, by conductivity, at large scales. 

Also magnetic flelds are affected by conductivity but to a lesser extent and only at sufficiently short 
scales (which are already erased by the flnite value of the thermal diffusivity scale, i.e. Silk damping). 
The typical magnetic diffusivity scale (i.e. the length-scale below which the magnetic fleld is dissipated 
by the flnite value of the conductivity) can be understood from the corresponding magnetic diffusivity 
equation. Inserting Eq. (12.461) into the second relation of Eq. (12.440 the magnetic diffusivity equation 
can be written as: 

^ = V X (4 X 5) + (2.47) 
or Ana 

According to Eq. (12.470 the magnetic fleld power spectrum will be diffused for typical wave-numbers 

k > kcr — \/Ano^H. Thus, only sufficiently short length-scales L < L„ c:^ are dissipated. The ratio 

of L„ to the Hubble radius, i.e. L^H, being suppressed by (T/Mp), is always minute. Around equahty 

we can estimate that L^jH ~ 3.9 x 10^^''(T/eV)^/^ where cTc is given by Eq. (I2.27p . 

For T < MeV the kinetic Reynolds number is smaller than one. This property is not verified, 
for instance, during the life of spiral galaxies where, effectively, the kinetic energy of the plasma can 
be converted into magnetic energy by means of the first term of Eq. (I2.47P which is often dubbed 
dynamo term [9|. When the kinetic Reynolds number is small (i.e. in the absence of kinetic turbulence) 
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the plasma description following from MHD can be also phrased in terms of the conservation of two 
interesting quantities, i.e. the magnetic flux and the magnetic helicity [7]: 



dt\Js 
dt\Jv 



Alia 



V X V X ■ dS, 



d^xB -V X B. 



(2.48) 
(2.49) 



Ana Jv 

In Eq. fl2.48p . S is an arbitrary finite surface that moves with the plasma. In the ideal MHD limit 
(i.e. a = aac oo) the magnetic flux is conserved. In the same limit also the magnetic helicity is 
conserved. In the resistive limit the magnetic flux and helicity are dissipated with a rate proportional 
to l/o" which is small provided the conductivity is sufficiently high. The term appearing at the right 
hand side of Eq. fl2.49p is called magnetic gyrotropy. Since, at high temperatures, the conductivity 
grows with T the ideal limit is always verified better and better as we go back in time. 

The conservation of the magnetic helicity is a statement on the conservation of the topological 
properties of the magnetic flux lines. If the magnetic field is completely stochastic, the magnetic flux 
lines will be closed loops evolving independently in the plasma and the helicity will vanish. There 
could be, however, more complicated topological situations where a single magnetic loop is twisted 
(like some kind of Mobius stripe) or the case where the magnetic loops are connected like the rings 
of a chain. In both cases the magnetic helicity will not be zero since it measures, essentially, the 
number of links and twists in the magnetic flux lines. The magnetic helicity will have no impact on 
our considerations since the scalar fluctuations of the geometry are not affected by the helical features 
of the magnetic fields. On the contrary in the vector and tensor cases the situation can be different 



In the resistive MHD approximation the electric components of the energy-momentum tensor can 
be neglected, while the magnetic components are present only at sufficiently large scales L > L^: 

52 . 1 



%\x,t) 



V(x.t) 



d2 

B,B^ - —6{ 
2 ' 



(2.50) 



where B^ = BiB^. In Eq. f l2.50p the contribution of the electric terms can be neglected since they are 
all suppressed by two powers of the conductivity. The Poynting vector can be also neglected since, at 
finite conductivity is suppressed as 



1 



-Ex B 



1 



-J X B, 



. ^ •■ (2.51) 

where Eq. f l2.44p has been used in the second equality. The spatial components of the energy momen- 
tum tensor can be phrased in terms of the magnetic pressure and of the anisotropic stress, i.e. 



%\x, r) = -6pb{x, t)61 + fll{x, r) 
where, with standard notations: 

B^(x) 



6pb{x, t). 



Spb{x, t) 



HKx,r) 



5pB 
3 ' 



B,B^ 



-Si 

3 ^ 



(2.52) 

(2.53) 
(2.54) 



It is practical to refer the magnetic fields to the photon background by means of the following rescaling 

6pb{x,t) _ B^{x,t) 



^b(x, t) 



Pyi 



T 



irrp 



d.dmi = {p, + p,)V'aB. 



(2.55) 
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where = a^p^. With the notations of Eq. fl2.55p the identity 

VVb = -^d,B,d^B'^ - iv^^B (2.56) 
allows to express the three-divergence of the Lorentz force appearing in Eq. (12.381) in terms of VL-q and 

4 4 

Stochastically distributed large-scale magnetic fields do not break the spatial isotropy of the back- 
ground geometry introduced in Eq. (12.11) . Nearly all magnetogenesis mechanisms suggest indeed that 
the large-scale magnetic fields should be stochastically distributed and characterized by their two-point 
function. Defining the Fourier amplitude of the magnetic fields as 

B,{x) = 777^ / dPkB,{k)e-'^-^, (2.58) 



(27r)3/2 

their two-point function can be expressed as 

{B,{k)Bj{p}) = '^P,^(k)PB{k)6^'\k + p), (2.59) 

where 

P^^ik) = (5,, - PB(fc) = ^b(|-)"'~'. (2.60) 

In Eq. (I2.60p ub is the magnetic spectral index and is the amplitude of the magnetic power 
spectrum referred to the magnetic pivot scale kj^. In terms of this definition the two-point function of 
the magnetic fields in real space can be written as 

dlnkP,,{k)PB{k)^^, r=\x-y\ (2.61) 

Different conventions exist in the literature for assigning the magnetic power spectrum. For instance 
in [251 the k~^ (appearing at the right hand side of Eq. fl2.59p ) was included in the definition of 
P-B{k). Those authors, indeed, only dealt with tensor and vector modes and were not confronted with 
the necessity of assigning power spectra according to the standards of CMB physics. We are forced 
on Eqs. (I2.58p - (l2.60p since these are the conventions used to define the power spectrum of curvature 
perturbations (see for instance [HI [T21 [HI [Il|) and it would be strange to use normalizations and 
definitions of the spectral indices that may differ from the ones which are commonly established when 
presenting theoretical and observational studies of the parameter space of CMB anisotropies. 



2.4 Gravitating magnetic fields 

All the species introduced so far gravitate and, therefore, affect the evolution of the metric perturba- 
tions and h. In Fourier space, the Hamiltonian and the momentum constraints stemming from the 

^Note that, according to Eq. (|2.58p . the dimensions of Bi{x) are of L^^. Thus, dimensionally, [Bi(k)] ~ L. But then 
it is easy to see (taking into account the dimensions of the three-dimensional Dirac delta function) that, dimensionally, 
[^b] = L~^, i.e. (and hence -Pb(^)) has the same dimensions of the magnetic energy density in real space. This is 
another good reason, unlike previous studies (see last three references), to follow the conventions expressed by Eq. 
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(00) and (Oi) components of the perturbed Einstein equations are, respectively 

2k'^^-nh' = 87TGa'^[6pt + 6pB], (2.62) 
= -'^7rGa^iPt + Pt)et. (2.63) 

In Eqs. (12.621) and (12.631) 6pt and 6t are the global density fluctuation of the plasma and the total 
velocity field defined as 

6pt = Spc + Sp^ + 6p^ + 6pb, (2.64) 

4 4 

(Pt + Pt)Ot = Y.(P^ + Pa)6'a = l^PuOu + l^P-^/O-/ + PcOc + PhOh- (2.65) 
a 6 6 

The spatial components of the perturbed Einstein equations (i.e., respectively, {i = j) and {i ^ j)) 
lead instead to: 

h" + 2nh' - 2k'^i = 2A7rGa'^[Spt + Spb], (2.66) 
{h + 60" + 2H{h + 60' - 2^^ = 2A7rGa^[{p, + p,)a, + {p^ + p^)(Tb]. (2.67) 

In Eq. (12.671) the neutrino anisotropic stress (also appearing in Eqs. (12. lip and (I2.12p ) has been 
consistently included. For analytical estimates it is also useful to write the evolution equation for the 
total density contrast which reads, in the synchronous gauge, 

Sp[ + m{cl + 1)8 p, + 37^5pnad + (Pt + Pt)Ot - {Pt + Pt)y = 0, (2.68) 

where the total sound speed c^^ and the non-adiabatic pressure fluctuation 5pnad have been introduced. 
Their respective definitions can be extracted from the following pair of relations: 

(5pt = 45pt + 5Pnad, 4 = 4- (2-69) 

Pt 

Equation (I2.69P implies that the pressure fluctuations can be generated either by inhomogeneities in 
the energy density or by fluctuations of the sound speed itself. The latter fluctuations are non-adiabatic 
in nature since they arise, physically, as a fluctuation of the specific entropy, i.e. the entropy density 
of the photon gas measured in units of the concentration of another given species. This property is 
customarily used to classify the initial conditions of CMB anisotropics which are therefore divided into 
adiabatic and non-adiabatic. 

In the case of the adiabatic mode, by definition, the fluctuations in the entropy density vanish over 
typical scales larger than the Hubble radius at recombination. The opposite holds for the non-adiabatic 
modes. In the case of the CDM-radiation mode the specific entropy is just given by <^ = T^/ric where 
Tie is the concentration of the CDM particles. The entropy fluctuations (i.e. the relative fluctuations 
in the specific entropy) are given by 

6^ 3 

More generally, given two species of the plasma the entropy fluctuations are defined as jl5l H6| W7\ 



S = — = -6^- 4- (2.70) 
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The non-adiabatic pressure fluctuations can then be written a^: 

5Pnad = pi Pj(Csi - cl^)Sij, c^^ = ^, (2.72) 

O'l-Pt ij Pi 

where Cgj and c^j are sound speeds of two (generic) species. Thus, according to Eq. (12.721) . 5pnad 
measures, indeed, the degree of compositeness of the plasma: if more species are present, more non- 
adiabatic modes are possible and 5pnad receives more contributions. In the pre-equality plasma there 
are four (regular) non-adiabatic modes (i.e. the CDM-radiation mode, the baryon-radiation mode, the 
neutrino density and the neutrino isocurvature modes). There is one adiabatic mode whose presence 
is strongly suggested by the analysis of cosmological data in the framework of a ACDM scenario. 
This does not exclude the presence of a dominant adiabatic mode with a subdominant non-adiabatic 
component (see, for instance, [5Tj). 



The analysis will be here limited to the case of a single adiabatic mode in the presence of large- 
scale magnetic fields. This is the minimal situation compatible with the ACDM framework. Our 
code can accommodate also non-adiabatic initial conditions in the presence of magnetic fields. The 
initial conditions discussed in [32], if appropriately translated to the synchronous frame, allow for this 
possibility. Mixed initial conditions involve a dominant magnetized adiabatic mode and a number of 
magnetized non-adiabatic modes with subdominant amplitude. 

Two variables are customarily used to parametrize the power spectrum of the metric fluctuations. 
They will be denoted by TZ and by C- In terms of the synchronous degrees of freedom, they can be 
defined as 

ri - ri Pt 
Even if both TZ and C, are gauge- invariant, their physical interpretation is obtained by expressing the 
two variables in specific gauges: TZ is often dubbed curvature perturbation since it corresponds, in 
the comoving orthogonal gauge, to the perturbations of the spatial curvature. In analog terms, C is 
interpreted as the curvature perturbation in the gauge where the density contrast vanishes (also called 
uniform density gauge). Taking the difference of TZ and C, and using the Hamiltonian constraint fl2.62p . 
the following equation can be obtained: 

2Pe -{h + Qj)' 

The quantity at the right hand side of Eq. (12.741) is 0{k!^T'^) and, therefore, it is negligible when the 
relevant wavelengths are larger than the Hubble radius, in particular around equality. This property 
can be immediately understood by expressing the combination at the right hand side of Eq. (12.741) 
in terms of longitudinal gauge variables and, most notably, if) which denotes the spatial fluctuation of 
the metric in the longitudinal gauge: 

In the longitudinal gauge is constant, to lowest order, when the relevant wavelengths are larger 
than the Hubble radius (see Appendix A). The evolution of ( on scales larger than the Hubble radius 
translates immediately in the evolution of TZ. The evolution of ( can be simply obtained by inserting 



^ Owing to their definitions, both Sij and Spnad are gauge- invariant. 
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Eq. (12.731) into Eq. (12.681) . The logic is to trade dpt in favor of (. The result of this manipulation, 
after the use of the covariant conservation equation of the total fluid (i.e. Eq. (I2.4p ). is: 



n 



pt + pt 



5Pnad + 



n 



Pt + pt 



h' + 6^' 



(2.76) 



where 6^ is the three-divergence of the velocity field in the longitudinal gauge (see the Appendix, 
Eq. (1A.2P ). If (5pnad = (as contemplated in the present paper), then Eq. (12.761) can be explicitly 
integrated with the result that 



(2.77) 



4(3a + 4) ' 

where a = a/aeq and = 1 — R^, (see Eq. (12.161) ). So the spectrum of the primordial adiabatic mode 
will be given in terms of C*(^) = 'Ti*{k). The two-point function in Fourier space will be, for TZ 



27r 



(7^,(A;,r)7^.(p,r)) = —Vn{k)5^^\k + p), Vn{k) = A 



(2.78) 



where rig is the scalar (adiabatic) spectral index, fcp = 0.002 Mpc~^ is the so-called pivot scale and An 
is, by definition, the amplitude of the power spectrum at the pivot scale. With these conventions the 
two-point function in real space becomes 

sin kr , ^ ^, 



/c 
d\nkVn{k)'- 



kr 



\x 



y\ 



(2.79) 



So far the evolution equations of the lowest multipoles of the Boltzmann hierarchy have been 
introduced. It is relevant to recall also the brightness perturbations of the radiation field which are 
related to the inhomogeneities of the Stokes parameters. In the synchronous coordinate system the 
evolution equations of the brightness perturbations can be written as 



A[ + ikfiAi = 

A'q + tkfiAQ 
A'u + ikfiAu 



-Ai + Aio + pv^ - ^P2{p)Sq 



-Aq + -(1-P2(/^))5q 



-e'A 



u, 



< + Hvk + -^{31 An + v^) + ,fc ^B 4aB 

-Kb 4/£b 



0, 



(2.80) 

(2.81) 
(2.82) 
(2.83) 



where R^ has been defined in Eq. (12.221) and where we defined Vh = 9\,/[ik). Moreover, in Eqs. (12.801) 
and (ESD): 

5q = Ai2 + Aqo + Aq2. (2.84) 

The notations A^ and Aq£ denote the £-th multipole of Aj and Aq. In Eqs. (I2.80p and (I2.8ip 
P2{.p) = (S/U^ — 1)/2 is the second Legendre polynomial. Equations (I2.80p - (l2.83p constitute the basis of 



the semi-analytical approach used to estimate the magnetized temperature autocorrelations [311 [32l [33] . 
In particular, the aforementioned equations have been solved in the tight-coupling approximation to 
first and second order [311 El]- The physical information contained in Eqs. fl2.80p - (l2.83p can be 
summarized by noticing that to zeroth-order in the tight-coupling expansion the CMB is not polarized 
in the baryon rest frame so that Aq and Au will be zero. To first-order in the tight-coupling expansion 
the quadrupole of the polarization (i.e. Aq2) is proportional to the zeroth-order dipole. Since the 
zeroth-order dipole feels the Lorentz force, the polarization is also affected by the presence of large- 
scale magnetic fields. 
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3 Magnetized initial conditions for the Boltzmann hierarchy 



At early times, close to the moment when initial conditions are set, the evolution equations for baryons 
and photons are integrated in the tight-coupling approximation. Otherwise this would represent a stiff 
problem owing to the largeness of the Thompson rate. Consider, first of all, the difference between 
the baryon velocity equation (i.e. Eq. (12.381) ) and the photon velocity equation (i.e. Eq. (I2.4ip ): the 
result of this manipulation is: 



Q^)' + — (1 + i?b)(^^7 - ^b) = fc'^ + ne,, - -— (i^B - 4aB), (3.i: 

itb 4 4itb 

where Eq. (I2.57P has been also used to express the Lorentz force in terms of fie and ctb- Owing to the 
presence of e', Eq. (13.11) stipulates that any initial difference in the baryon-photon velocity is quickly 
washed out. Consequently, at early times 6*^ ~ ^^b- 

Denoting by the common value of the photon-baryon velocity field, the corresponding evolution 
equation can be obtained by combining Eqs. (I2.38P and (I2.4ip in such a way that the scattering terms 
exactly cancel at the price of introducing explicitly i?b i-e. the baryon-to-photon ratio of Eq. (12.22p . 
The net result of this procedure is: 

' 1 + i?b p^{Rh + 1) 4(1 + i?b) 4(1 + i?b) 

where we have also taken into account the shear viscosity contribution (proportional to rf) which is 
responsible of the diffusion damping: 

4 1 
V = —P'y>^Th, Axh = -• (3.3) 

15 e' 

The shear viscosity term (to a given order in the tight-coupling expansion) allows for the estimate of 

diffusive effects. Standard considerations related to the zeroth-order in the tight coupling expansion 

imply that 

1 2 r . aodr' 



,2 - .1 Csb(r')^V^ . (3.4) 

To second order in the tight-coupling expansion the inclusion of the polarization allows to estimate 



dr' 



16 ^ B 



15 i?b + 1 



(3.5) 



kl Jo 6(i?b + l)e' 

The factor 16/15 arises since the polarization fluctuations are taken consistently into account in the 
derivation. This difference is physically relevant. Grossly speaking we can indeed say that more 
polarization implies more anisotropy (and vice versa); more polarization implies a faster damping by 
diffusion. Note that k-o provides an effective ultra-violet cut-off for the magnetic energy spectra and 
will be used later on. 



= - -6^1,, ^h = — - ^7b- (3.6) 



Correspondingly Eqs. (12.370 and (I2.4ip can be written as 

2 4 h' 
-h' - -9,^, 61, = - 

Equations (13. 2p and (13.60 can be further combined to get a single equation for the density contrast of 
the radiation field 6^ with the result that 

HR^ 2 2 e 2 [CRb + lM _^ o 1 (-^7^ 
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where the Silk damping has been neglected. 

The same equation can be obtained by systematically expanding the brightness perturbations of 
the radiation field (see, in particular, Eqs. (12.801) and (I2.83P ) to zeroth order in the tight coupling 
expansion and by recalling that the precise relation between the monopole of the radiation intensity 
and the photon density contrast is given by 4Aio = 5^. In Eq. (13.71) the baryon-photon sound speed 
Csb has been introduced (see Eq. (12.211) ). In the absence of magnetic fields, the second source term 
in Eq. (13.71) vanishes. The resulting equation (in different gauges) can be again employed for the 
semi-analytical estimates of the temperature autocorrelations [531 EH EH EH EZ] • 

The whole Einstein-Boltzmann hierarchy will now be solved to zeroth-order in the tight-coupling 
expansion and for typical wavelengths larger than the Hubble radius before equality (i.e. kr < 1 for 
T < Teq)- The Hamiltoniau constraint of Eq. (I2.62p and Eq. (I2.66P will first be solved. The obtained 
solution, parametrized in terms of a suitable number of arbitrary constants, will be inserted into the 
other equations. The final solution will only depend upon the spectrum of the adiabatic mode and 
upon the power spectra of and ctb- Equations (I2.62p and (I2.66P are solved provided: 

^{k, t) = -2C{k) + A^{k)k^r\ h{k, r) = -C{k)eT^ - Ah{k)k^T\ (3.8) 
S^{k, r) = -R^riBik) - A^{k)eT\ 6^{k, r) = -R^QB{k) - A^{k)k^T^. (3.9) 

To lowest order in kr, h{k, r) does not have a constant term whose presence would entail a spuri- 
ous gauge mode which must be projected out by exploiting the remaining gauge freedom [38]. The 
compatibility of Eqs. (13.80 and (13. 9p with Eqs. (12.620 and (I2.66P leads to the following condition 

3{R^A^{k) + RuA^{k)) = 2C{k), (3.10) 

which guarantees that Eqs. (I2.62p and (I2.66P are satisfied with corrections which are 0{k'^T'^). In 
Eq. ( ]3.10p R^ = 1 — Ru (see Eq. ( 12.160 ). The neutrino and photon fractions arise since, prior to 
equality, the Hamiltonian constraint of Eq. (12.621) can be simply written, after using Eq. (12.21) . as 
4A;^^ — 2T-Ch' = 3T-C^[R^6u + R^6^]. Equation (I2.66P can be also recast in a similar form. 

The evolution equations of the velocity fields can be solved with a similar technique. In particular 
Eqs. (1230D - (|2TTD and ([S^D-dMD imply that 

e^^{k,T) = D^i{k)k\ + D^^i{k)k''T\ e,{k,T) = D^^\k)k\ + D^;^\k)k'T\ (3.11) 

where the various constants must satisfy: 

3A^(fc) = 2D^:^i{k) + 2C(fc), 3A,(fc) = 2D\^\k) + 2C(fc), (3.12) 

<b^(^) = X^b(^) - ^^^^)' = - ^^^(k). (3.13) 

The compatibility of the obtained solution with the momentum constraint of Eq. (12.630 imposes, 
moreover, the following pair of conditions 

RuD^J\k) + R^D^^iik) = 0, RuDf\k) + R-yD^^{k) + A^{k) = 0. (3.14) 

The evolution equations involving the neutrino anisotropic stress (i.e. Eqs. (I2.12p and (12.670 ) are 
solved by 

(^Ak, r) = -^Mk) + A„{k)eT\ (3.15) 
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As in the previous cases there are non-trivial conditions to be satisfied and they are, in the case of 
Eqs. (12121) and ([2l7D, 

3[6A^{k) - Cik)] + 2C{k) = GKA^ik), A^ik) = ^D^^\k) + ^C{k) - \A^{k). (3.16) 

15 15 5 

We are left with the evolution equations of the baryon and CDM density contrasts and with the velocity 
field of the CDM perturbations. The solutions for these quantities, as they emerge, respectively, from 
Eqs. dMD and ([M]) are: 

■\R,n^{k) - ^fcV, 

3 „ ^ 1 r^.,^ R. 



ec{k,T) 



0, 



Sc{k,T) 



6b{k,T) 



--R,n^{k) - - 



C{k) - a^ik) + ^n^ik) 



k^T^ 



(3.17) 



All the compatibility conditions constraining the form of the solution can be solved, and, after some 
algebra, the full solution for the initial conditions of the lowest multipoles of the Einstein- Boltzmann 
hierarchy becomes: 



h{k, T 



Sj{k, T 
5u{k,T 

Sc{k,r 
Sb{k,T 
9yhik,T 

9,{k,T 

ec{k,T 
cr^{k,T 



-2C{k) + 



4R, + 5 ^^j^^^R.,{AaB{k)-RMk)) 



-C{k)k'r' - 



6(4i?, + 15) 
1 



6(4i?, + 15) 



2_2 



k'r 



(3.18) 
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36 
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36 
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k^T\ 



(3.19) 

(3.20) 
(3.21) 

(3.22) 
(3.23) 
(3.24) 

(3.25) 
(3.26) 
(3.27) 



3(4i?^ + 15) 2R^{AR^ + 15) 

In the limit flB^k) — * and (Tb(A;) — > the solution corresponds to the usual adiabatic mode. To 
lowest order, in fact, we can appreciate that the non-adiabatic pressure fiuctuations introduced in Eqs. 
( ETTj) and fl^T2]) vanish since 



4 4 
6^{k,T) ~ S^{k,T) ~ -5c(A;,r) ~ -6b{k,T). 



(3.28) 



The adiabatic mode is parametrized in terms of C{k). The spectrum of C{k) is related to the spectrum 
of 7?-*(/c) introduced in Eq. (12.781) . Equation fl3.18p can be inserted into Eq. (12.731) with the result that. 
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to lowest order in kr, lZ{k) = —2C{k). According to Eq. (13.261) . the CDM velocity field is set exactly 
to 0. This requirement avoids the presence of the second gauge mode (the first one was projected 
out by excluding the constant solution for h{k,T)). For infinitesimal diffeomeorphisms (see also the 
Appendix) the metric fluctuations change by the Lie derivative in the direction parametrized by the 
two gauge parameters eo and e. By imposing the synchronous gauge the eo and e are determined only 
up to two functions which are constant in time but not in space. This is the ultimate rationale for 
the possible persistence, in the synchronous gauge, of two spurious gauge modes. Such a possibility 
is avoided by gauging away the two unphysical solutions and this achieved, at a practical level, by 
setting 6c{k,T) = and by requiring that h{k,T) does not have a constant mode. A complementary 
way of addressing this issue is to work out the same solution in a frame where the gauge freedom 
is completely removed. The solution obtained in this gauge must match with the results obtained 
in the synchronous gauge and transformed to the new gauge. In the case of the present problem, 
also following previous works, it is productive to cross-check the results in the longitudinal coordinate 
system. The main ingredients for this analysis are reported in the Appendix where the synchronous 
and longitudinal solutions are explicitly connected. 

The initial spectrum of the magnetic fields is encoded in and (Tb. The spectrum of Qsik) and 
crB(^) must then be computed in terms of the spectrum of the magnetic field introduced in Eqs. (12.591) 
and (I2.60p . This is a rather lengthy calculation and here only the main steps will be outlined. 

Since both Qb{x) and o"B(a;) are quadratic in the magnetic field intensities (see, for instance, Eqs. 
(I2.55P and (I2.56P ) their expressions in Fourier space will lead to the following two convolutions 

^ ^ '' d^kBi{k)B\q-k), 



(27r)3/2 Svrp, 



(27r)3/2 IQTip 



(fk 



3{q^ - y)k' 



B,{k)Bi{q-k)-Bi{q-k)B\k) 



(3.29) 
(3.30) 



The correlation functions for VLB{k) and (JB{k) are then defined as 



27r 



27r2 



(3.31) 



To compute Vn{q) and Va{q) in terms of the magnetic power spectra we must go through the straight- 
forward (but rather lengthy) procedure of expressing the stochastic averages of four fields in terms of 
the two-point function of Eqs. (12.591) and (I2.60p . Then the obtained results must be integrated over 
the momenta. After performing the first of the previously mentioned steps we obtain 
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(3.32) 



1 + 



k'^\p — /cp 



k'^\p- 
2[k- (p 



k)Y- 



(3.33) 



In Eqs. (I3.32P and (13.330 the notation k = \k\ has been employed. To complete the calculation the 
angular integration and then the radial integration must be performed. The integration measure can 
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be written, in spherical coordinates, as d^k = k'^dkd cosddip. The integration over cp is trivial and leads 
just to a factor 27r. However the integration over dcosd (between —1 and 1) is rather cumbersome 
(but doable in exact terms). Indeed all scalar products arising in Eqs. ( 13.32^ and (13.330 induce a 

— * 

factor cos'd. So calling x = cost?, the expressions like \q — k\ (which appear ubiquitously in Eqs. (13.321) 
and (13.331) ) become \q — k\ = ^Jcp- + k'^ — 2qkx. The combinations of all similar factors must then be 
integrated over x. The physically interesting region of spectra is realized when 1 < < 5/2. We 
will conventionally refer to this case as to the one of blue spectral indices. It is also interesting to 
discuss in some detail the case of red spectra (i.e. ne < 1). Violet spectra (i.e. ti-q ^ 1) are mainly 
constrained by the diffusion scale kj) (see Eq. (13.51) ). Simplistic estimates of the Silk damping scale 
lead to /c^^ — 0.3(a/adec)^''^/[^A^'^b] Mpc^. The exactly scale-invariant case leads to a logarithmically 
divergent power spectrum. 

The magnetic power spectra are usually defined within an appropriate regularization of the mag- 
netic energy density [30] : 



{B,{x)B\y)) = 2 / dlnkP^ikf-^Wik), 



\x - y\ 



(3.34) 



where W{k) is an appropriate window function. Consider first the case of a blue spectrum. The 
energy density can be regularized over a typical comoving scale L (which is related to the magnetic 
pivot scale in Fourier space) by means of a Gaussian window function W{k) = e-^ ^ . Equation fimi) 
then implies: 



Bl{r) = (27r)i-"MBr 



?2b — /riB — 1 3 r^/^L 



'11 



2 2 '2' levrV 

where Fii(a, 6, z) is the Kummer confluent hypergeometric function [l8lll9]. Since lim^^o -^11(0, z 
1, 



(3.35) 



El = \immx)Bj{y)) = AU27ry-^-Ti 



(3.36) 



In the radial integrals of Eqs. (13.321) and (13.331) . ^b can be traded for so that VniJ^) and Vaik) 
can be written, respectively, as: 
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(3.37) 

(3.38) 
(3.39) 

(3.40) 



It should be remarked that when 1 < riB < 5/2, we can formally send the diffusion scale to infinity 
(i.e. /cd — > 00) and the final result will still be convergent. Consequently, as already remarked in 
related contexts [28] the diffusion damping only enters the case when the spectral slopes are violet (i.e. 
ub > 5/2). 
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For ne < 1 the window function appearing in Eq. fl3.34p can be chosen as a simple step function 
W{k) = 6{k — ko). If this is the case and Vaik) can be formally written exactly as in Eq. (13.371) 

but with two slightly different pre-factors which shall be denoted by T{n-Q) and Q{n-Q): 
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(3.41) 



nB) 



(riB - l)(2nB - 5) 1 - ub 
Aub + 116 8 1 



15(5 — 2nB)(nB 



+ 



( kn 



nB-l' 



ne-l 



3 1 — TiB V ^ 



(3.42) 
(3.43) 



where k^ is of the order of (but smaller than) the Hubble rate. 

For normalization purposes it is useful to have an explicit expression of the Sachs- Wolfe plateau 
which includes the magnetic energy density. This estimate can be performed by solving Eq. (12.801) 
with the line integration method. The result of this procedure is 
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where /C(r) is the visibility function and 
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(3.44) 



(3.45) 



The term jj? appearing in Eq. (13.441) can be integrated by parts and, subsequently, the visibility 
function can be approximated by a Dirac delta function centered at the decoupling time. Neglecting 
the integrated Sachs- Wolfe contribution and the Doppler term: 



,(SW) 
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5, jh + GQ" 
4 2P 



(3.46) 



Tdo 



where T^ec has been neglected in comparison with tq in the argument of the exponential factor. To 
evaluate Eq. (13.461) we need to know the value of the combination {h + 6.^)' after equality when the 
relevant modes have wavelengths larger than the Hubble radius. Let us notice that, for the mentioned 
wavelengths, Eq. (13.61) implies that S'^ ~ 2/i'/3. Thus, denoting by S!:^^ and 5^*-' the final (i.e. at the 
decoupling) and initial (i.e. before equality) values of the density contrast we will have 



6? = 6^^ + ^{h(^^-h(^^). 



(3.47) 



The evolution of ^ across equality can be obtained from Eqs. (I2.73p . (12.761) and (12.771) and it is given 
by solving the following equation: 
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(3.48) 
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where a = a/a^q. Once ^ is known we can easily deduce {h + 6^)' from Eq. (12.671) which imphes 
(neglecting the anisotropic stress when the corresponding wavelengths are larger than the Hubble 
radius) : 

[{h + 60V]' ~ 2ea^^. (3.49) 
The final result for the ordinary Sachs- Wolfe term can then be written as 



n^{k) R 



+ 



-ikjjiTQ 



(3.50) 
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Expanding the plane wave in series of Legendre polynomials the ^u{k^ tq) can be easily extracted and 
the estimated with standard integration over the comoving wave-number k. The result is: 
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(3.54) 



In Eq. fl3.54p 7br is the correlation angle that has been included to keep the expressions as general 
as possible. In what follows the main focus will however be on the case where the adiabatic mode of 
curvature perturbations is not correlated with the magnetized contribution (i.e. 7br = 7r/2). Note, 
however, that if cos7br > then the cross-correlation between the adiabatic component and the 
magnetic component will lower the Sachs- Wolfe plateau allowing for a magnetized contribution which 
is comparatively larger than in the case where 7br = vr/2. We leave this possibility for future studies 

Em. 



4 Magnetized temperature autocorrelations 

In Fig. [1] the results of the numerical integration are illustrated in terms of the temperature autocor- 
relations. The parameters are fixed to the best fit of the WMAP data alone implying that the value 
of the scalar spectral index is ris = 0.958. The full set of cosmological parameters used to compute the 
models in Fig. [T]is given as follow^: 

(l^bo, ^co, ^A, ^0, ^^s, r) = (0.042, 0.198, 0.76, 0.732, 0.958, 0.089). (4.1) 

^Consistently with our notations we should denote with e^e the optical depth to reionization which is conventionally 
denoted by r and which we use to indicate the conformal time. However, since in the present and in the following sections 
the optical depth and the conformal time coordinate are never mentioned together, we will stick to the conventional 
terminology and denote with t the optical depth. 
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Figure 1: The temperature autocorrelations for blue magnetic spectral indices are compared with the 
best fit model arising from the WMAP alone analyzed in terms of a pure ACDM model with no tensors. 
The value of the magnetic pivot scale is /cl = 1 Mpc~^. 

Moreover the tensors are absent from the fit and Ri, = 0.408 (i.e., according to Eq. fl2.16p . N^, = 3.04). 
In Fig. [1] (full curve in both plots) the Cj,^"^^ are illustrated for the best fit parameters reported in Eq. 
(14. ip . The magnetic spectral indices are, in both cases blue, i.e., according to the terminology of the 
previous section, 1 < n-Q < 5/2. With the dashed line the magnetic fields corresponding to i?L = 0.1 
nG is reported. The dashed curve cannot be distinguished from the best fit curve. If the regularized 
magnetic field intensity is i?L < (9(0.1 nG), then the difference of the TT correlations (as well as the 
EE and TE correlations) with respect to the three year best fit is below the accuracy of the code. The 
latter statement depends, of course, on the spectral index and on the range of multipoles. Indeed, as 
argued in Section El the large-multipole region (i.e. £ ^ 1500) is more sensitive to regularized fields 
of nG strengths. 

The dot-dashed curve denotes, in both plots, the temperature autocorrelations computed in the 
case -Bl = 50 nG. The results illustrated in Fig. [1] are qualitatively similar for different choices of 
the parameters close to the best fit values. The inclusion of a magnetized background has a threefold 
effect on the temperature autocorrelations. The height of the first peak gets increased. The second 
peak is distorted and it eventually turns into a hump for sufficiently large values of i?L (or of ne). The 
third peak is, at the same time, distorted and raised. In Fig. [T] as we move from the plot at the left 
to the plot at the right the spectral index increases. The increase of the spectral slope entails also an 
increase in the distortions. The latter trend, however, is not monotonic at least in the case of the first 
acoustic peak. A more thorough illustration of this feature will be provided in Section [61 It should 
be borne in mind that, within the conventions established in Section [3] the scale-invariant limit of the 
magnetic power spectra is realized for ne — ^ 1. In analog terms, the Harrison- Zeldovich limit for the 
power spectrum of curvature perturbations occurs when ng —>■ 1. 

The features illustrated in Fig. [1] and scrutinized in the previous paragraph do not depend upon 
the data sets. The same qualitative patterns can be observed if the pivotal model is taken to be the 
best fit inferred from the combination of the WMAP data with all the other data. In this case the 
central values of the cosmological parameters are slightly changed [HI |T2l US] according to: 

(n^o, n^o, ^A, ho, n„ r) = (0.044, 0.223, 0.733, 0.704, 0.947, 0.073). (4.2) 
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Figure 2: The temperature autocorrelations for blue magnetic spectral indices. In the plot at the left 
the cosmo logical parameters are fixed to the central values of the best fit when the WMAP data are 
combined with all the cosmological data sets (see Eq. fl4.2l) ). In the plot at the right the cosmological 
parameters are fixed to the central values of the best fit when the WMAP data are combined with the 
gold sample of type la supernovae (see Eq. (14. 3p ). As in Fig. [T]the pivotal model is ACDM scenario 
with no tensors. The value of the magnetic pivot scale is /cl = 1 Mpc~^. 

By fitting the WMAP data with the ones of the gold sample of type la supernovae [20] the central 
values of the cosmological parameters are yet a bit different from the ones reported in Eqs. (14.11) and 

(fibo, ^co, ^A, ho, n„ t) = (0.045, 0.231, 0.724, 0.701, 0.946, 0.079). (4.3) 

In Fig. [2] the temperature autocorrelations are computed when the cosmological parameters are fixed 
as in Eq. (14. 2p . As the magnetic field strength increases from 10 to 50 nG the distortion patterns 
already illustrated in Fig. [1] become more pronounced. By comparing Figs. [1] and [2] the same distortion 
patterns can be observed. In Fig. [2] the spectral tilt increases from the left to the right plot. As in 
Fig. [H also in Fig. [2] the plot at the left is close to the scale-invariant limit of the magnetic power 
spectrum while the plot at the right illustrates a spectral slope which is bluer. 

The distortion patterns arising in Figs. [1] and [2] have a semi-analytical interpretation. In [33] the 
effects of magnetic fields on the temperature autocorrelations have been discussed in a semi-analytical 
perspective and for blue spectral indices (i.e. 1 < riB < 5/2). A consistent use of the tight-coupling 
approximation allowed for the estimate of the at low multipoles (i.e. i < 30) and also at large 

multipoles (i.e. i > 100). Using the large-order expansion of the spherical Bessel functions (and of 
their derivatives) the shape of the TT correlations has been reduced to the numerical calculation of 
four integrals [33] . The semi-analytical approach described in |33] seems, a posteriori, rather brutal. 
Nonetheless it is amusing that the essential patterns of the distortions in the acoustic region have been 
correctly captured. 

The Doppler region is sensitive to the relative phases and amplitudes of the monopole and dipole 
terms of the brightness perturbations. The contribution of the Lorentz force provides an extra source 
to the monopole equation (see Eq. (13.71) ). A computable difference in the relative amplitudes of the 
monopole and dipole terms then arises also analytically and it is reflected in the overall distortion. 
The reasonable results obtainable within the tight-coupling expansion (improved to second order) 
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Figure 3: The temperature autocorrelations are illustrated in the case of red tilt. The values of the 
cosmological parameters in the left and right plots have been fixed, respectively, as in Eqs. ( 14.2p and 

dH). 

represent a powerful cross-check for the consistency both of the numerical approach and of the semi- 
analytical calculation. In spite of the encouraging agreement of the numerical results with the semi- 
analytical evaluations, a weaker dependence upon the amplitude of the magnetic fields has to be 
admitted. Indeed, the reionization effects have been neglected in [33] and the recombination has been 
treated within a Gaussian parametrization of the visibility function. Furthermore, always in [SSj, the 
overall amplitude of the angular power spectrum was determined by matching the low-i regime with 
the Doppler regime where the Bessel functions have been basically replaced with their asymptotic 
expressions for i ^ 1. Numerically these approximations have been dropped. 

The range of spectral indices 1 < < 5/2, on a theoretical ground, is well motivated. The 
two-point function of the magnetic fields decreases, in this case at large distances. The diffusive 
effects are negligible since, as explained in the previous section, the two-point functions of the energy 
density (and of the Lorentz force) is insensitive, in this case, to the ultra-violet cut-off. Finally various 
magnetogenesis models predict this kind of spectral^. 

The last motivation, however, is just accidental. The true question behind these considerations is 
slightly different and can be phrased by asking: which is the spectrum of magnetic fields at the onset of 
gravitational collapse of the protogalaxy? As we can in principle measure the matter power spectrum 
it would not be insane to think that, in a future, also the magnetic power spectra of different objects 
could be measured. Indeed there are attempts to characterize, for instance, the present features of 
our galactic magnetic field in terms of an appropriate power spectrum [59]. In a related perspective 
one could observe that it is equally plausible to study the mean-squared fluctuation of the Faraday 
Rotation Measure (RM), as it was proposed in Refs. [60], [61]. One of the key projects of the radio- 
astronomy community is the celebrated Square Kilometer Array (SKA)0. This instrument will allow to 

^°See, for instance, |9j and [SO] for some reviews on this subject. It would be impossible to refer to all the attempts 
along this direction. For recent results see [58] and references therein. As specifically discussed in the introduction, 
the purpose her is not to endorse a particular model but to develop the tools which will allow to assess the primordial 
nature of the magnetic field. In this sense the goal of the present analysis is more modest. 

^^The collecting area of the instrument, as the name suggest, will be of 10^ m^. The specifications for the SKA 
require an angular resolution of 0.1 arcsec at 1.4 GHz, a frequency capability of 0.1-25 GHz, and a field of view of at 
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Figure 4: The temperature autocorrelations are illustrated in the case of vanishing adiabatic mode 
for fixed magnetic field strength (plot at the left) and for fixed magnetic spectral index (plot at the 
right). The magnetic pivot scale is taken to be /cl = lMpc~^ and the other values of the cosmological 
parameters are fixed to their best fit values as in Eq. f l4.2p (plot at the left) and as in Eq. ( 14. 3p (plot 
at the right). 

obtain full sky survey of the RM and, in that context, it will be even more plausible to collect valuable 
informations on the magnetic power spectra at large scales. These measurements, even if feasible in 
the future, will not provide direct indications on the protogalactic field but rather on the present field. 
Still it is not excluded that the morphological features of the observed field could be connected with the 
protogalactic features. If the magnetic field does not flip its sign from one spiral arm to the other, then 
a strong dynamo action can be suspected [61] . In the opposite case the magnetic field of spiral galaxies 
should be primordial i.e. present already at the onset of gravitational collapse. An excellent review on 
the evidence of magnetism in nearby galaxies can be found in [65j. In the model- independent approach 
followed in the present paper it is natural to ask what happens if the magnetic power spectra have a 
red tilt. In Fig. [3] (left plot) the spectral tilts of the magnetic power spectrum and of the spectrum of 
curvature perturbations coincide. In other words Ug = ub = 0.947 < 1. The specific figure (i.e. 0.947) 
is dictated by the adoption, as fiducial set of data, of the best fit to the WMAP data combined with 
all the other data (see Eq. (14.21) and also the titles of the plots in Fig. [3]). In the plot at the right 
of Fig. [3] the spectral indices of the magnetic energy density and of the anisotropic stress decreases 
(i.e. ub = 0.5 <C Us). So we can say that, in Fig. [3] the magnetic spectral index is redder at the right 
than at the left. The two plots in Fig. [3] have been presented in semi-logarithmic coordinates and 
this choice allows to scrutinize in more depth the main feature associated with red spectral indices: 
as the spectral index becomes redder, a systematic decrease of the low multipoles is observed. This 
trend has been investigated with other examples (which will not be reported for reasons of space). 
The conclusion of this analysis is that, indeed, red spectra can lead to a lower quadrupole and, only 
apparently, improve the agreement with the data. We say only apparently because every time the TT 
correlation diminishes at low multipoles, the first pair of acoustic peaks is raised and distorted to an 

least Ideg^ at 1.4 GHz [5^ [55], The number of independent beams is expected to be larger than 4 and the number 
of instantaneous pencil beams will be roughly 100 with a maximum primary beam separation of about 100 deg at low 
frequencies (becoming 1 deg at high frequencies, i.e. of the order of 1 GHz). These specifications will probably allow 
full sky surveys of Faraday Rotation. 
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Figure 5: The effect of tlie magnetic fields is illustrated in the case when the CDM fraction is increased 
by enforcing the spatial flatness and when the baryonic fraction is kept fixed to the best fit value. 

unacceptable degree. This aspect can be appreciated, in a rather extreme case, in the plot at the right 
of Fig. [3l For i?L = 50 nG the low multipoles would indeed represent better the experimental points. 
However, the Doppler peak explodes to, roughly, 9000 (/xK)^. 

When the spectra have a red tilt the magnetic pivot scale coincides effectively with the infra-red 
cut-off of the spectrum. If the cut-off is as large as the present Hubble patch the magnetic field 
acquires, for practical purposes a preferred direction. One could be tempted to say that this offers 
an explanation of the lower value of the quadrupole. Indeed various proposals have been put forward 
to explain the quadrupole with a specific anisotropic model falling in one of the Bianchi classes. The 
consistency of our numerical approach shows, however, that these kinds of red spectra lowering the 
quadrupole are simply pathological. Indeed, as it appears from Fig. [3l red spectral slopes lowering the 
quadrupole and spontaneously breaking spatial isotropy are not consistent at higher multipoles. We 
suspect that current attempts of justifying the quadrupole anomaly based on anisotropic models would 
simply fail when confronted with the higher multipoles. The considerations reported here are a first 
quantitative indication in this direction. It would be interesting to pursue this analysis in further detail 
[50] . The quadrupole anomaly is probably better addressed in the framework of pre- inflationary initial 
conditions such as the ones discussed in [BSl EZ] where the low value of the quadrupole is attributed 
to a fast-roll phase. 

The final comment in connection with red spectra deals with regularized magnetic field of regu- 
larized amplitude much smaller than the nG. In both plots of Fig. [3] the case of -Bl = 0.1 nG cannot 
be distinguished from the corresponding best fit model. The latter observation suggests that weak 
magnetic fields (i.e. -Bl — 0(0.1 nG)) with red tilt are not incompatible with current data on the TT 
correlations. 

If the adiabatic mode of primeval origin is switched to zero, we can expect, on the basis of purely 
analytical considerations, that the typical amplitude of the temperature autocorrelations will be van- 
ishingly small to begin with. Indeed, the overall amplitude will not be controlled by the power spectrum 
of the adiabatic mode but by the power spectrum of the magnetic energy density, i.e. Vn{k). The am- 
plitude of the adiabatic power spectrum evaluated at the pivot scale kp is of the order of 10"^. The am- 
plitude oiVn{k), evaluated at the magnetic pivot scale /cl goes as 

^BL which is C(lO-iS) for 5l ~ nG. 
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Consequently if the adiabatic contribution is switched off we will expect that the temperature auto- 
correlations will be about 9 orders of magnitude smaller than in the case where the adiabatic mode 
was present. This means that while in the case of, for instance, Fig. + 1)C^^'^'' / {27r) ~ 10^(/iK)^, 
when the adiabatic mode is switched to zero i{i + 1)C^^'^^ / (2^) ~ 10~^(/iK)^. 

In Fig. m the theoretical expectations are confirmed by the numerical results. In the plot at the 
right i?L = 1 nG. For different values of the spectral indices the TT correlation exhibit a humpy profile 
in the acoustic region. This point is also stressed in the right plot of Fig. H] where for fixed tt-b the 
magnetic field strength is enhanced from 1 nG to 5 nG. This entails an increase of the temperature 
autocorrelations of a factor 6 x 10^ which fits with our expectation which would be, in this case, 
(5)^. According to Fig. [T]the case -Bl = 50 nG is already excluded by the present data and for the 
corresponding values of the spectral index. This observation helps along two opposite directions. The 
nature of the distortion induced by the magnetic fields seems to be hard to reproduce by varying 
the standard CMB parameters. For instance it is known that by lowering cUb the height of the peaks 
diminishes. Similar effects (but with a different quantitative impact) are observed when Uc increases 
(always enforcing the flatness of the model). None of these two effects distorts the peaks as in the 
case of magnetic fields. One can also think that by adding spatial curvature and by either decreasing 
tUb or increasing the effects of the magnetic fields can be appropriately mimicked. However, in the 
case of the magnetic fields not only the first peak increases but also the ratio of the second to the first 
peak is modified. The shift in the position of the peaks is much more severe in non-flat models than 
in the case of nG magnetic fields. 

In Fig. [5] (plot at the left) the baryonic fraction f2bo has been fixed to the best fit value of the 
WMAP data alone (see Eq. fl4.ip ). The CDM contribution has then been increased (by always 
keeping the model flat). In the right plot we took the most extreme model illustrated in the left plot 
(i.e. the one labeled by the dot-dashed line) and compared it with the same model where, however the 
magnetic field is included. This shows that the kind of correlated distortion induced by the magnetic 
fields cannot be simply reduced to an increase of the peaks (see also Section E] for a more extended 
scrutiny of this statement). We analyzed the characteristic shapes obtainable by changing also other 
parameters as the Hubble rate and, also in that case, the magnetic fields induce distortions which 
cannot be mimicked by known shape effects. Under certain circumstances, a slight increase in the 
CDM fraction (less extreme than those illustrated in Fig. [5]) can be compensated by the presence of 
a minute magnetized background. 

This type of considerations bring up naturally the need of including the magnetic fields as an 
extra set of parameters in the current strategies of data analysis [25J. In its simplest realization the 
magnetized ACDM paradigm entails the inclusion of two new parameters, i.e. the magnetic spectral 
index and the amplitude of the regularized field. It cannot be excluded, in other words, that a combined 
action of different effects will be compensated by a magnetic field leading, ultimately, to a better fit. 
The accuracy of forthcoming data (see also Section [6]) seems to suggest that we will soon be sensitive 
to nG magnetic fields and, then, global strategies of parameter extraction will allow either to confirm 
the ACDM paradigm or to improve it. 

A handy parallel can be drawn with a slightly different physical case which however bears some 
analogy with the one discussed here. A commonly employed approach to the initial conditions is the 
one we could define as model- independent (see, for instance, [511 [68]). When analyzing cosmological 
data a very interesting question is to know if the data allow for a sizable non-adiabatic component. 
It is by now well established how to constrain the CDM-isocurvature mode [BSl EO]- This analysis 
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Figure 6: The magnetized TE correlations are illustrated for different values of the spectral indices 
and different values of the magnetic field intensities. The fiducial set of the cosmological parameters 
is the one of Eq. fl4.ip . 

entails, in the simplest case, the addition of two extra-parameters, i.e. the amplitude and spectral 
index of the non-adiabatic mode. However there could be even more complicated situations where a 
cross-corrrelation term is present. This term parametrizes a possible correlation between the adiabatic 
an the non-adiabatic components and typically leads to further parameters. Recently interesting 
results have been reported in this context. For instance it has been shown that indeed the addition of 
an adiabatic component with blue spectrum may improve the global fits of cosmological parameters 
[71] . The situation described in the case of the isocurvature modes is similar to what happens in the 
case of large-scale magnetic fields with the crucial difference that, in the present case, not only the 
initial conditions but also the dynamics is affected by the addition of stochastic magnetic fields. As 
customarily done, for other parameters, it will be appropriate to include the magnetic fields when 
confronting all the cosmological data sets. This idea will allow to set bounds and compare fits in a way 
which is less brutal than the one sometimes employed when dealing with large-scale magnetic fields. 
The numerical approach developed and applied in the present study is the first step in this direction 
which we plan to investigate throughly in the near future [5U] . 

Always in connection with the isocurvature modes we wish to stress that large-scale magnetic fields 
can be included also in the case when the initial conditions are not predominantly adiabatic but rather 
obtained as a mixture of adiabatic and non-adiabatic components. In this study, for reasons of space, 
we just focused on the magnetized adiabatic mode. It is therefore possible to study, with our approach, 
all the usual situations encountered in conventional CMB calculations |50j . 



5 Polarization correlations and cross-correlations 



The TE cross-correlations are probably the strongest indicator of the adiabatic nature of the CMB 
initial conditions. Indeed, in the adiabatic case, the Cf^'^^ shows a characteristic anticorrelation peak 
for i ~ (3/4)£Doppier — 150 where ^ooppier denotes the observed position of the Doppler peak [HI [32] . 
The relation between Doppler and anticorrelation peaks is a distinctive feature of the adiabaticity of 
the fluctuations prior to recombination. The polarization observations are therefore a rather sensitive 
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Figure 7: The magnetized EE correlations for blue values of the magnetic spectral tilt. The cosmo- 
logical parameters are the same as in Fig. O 

tool which can be used to scrutinize the possible contribution of a magnetized component. At the 
moment various experiments reported a positive detection of the EE and TE correlations. Besides 
the three year results of the WMAP collaboration [12] there are, at the moment, the three year 
data of the DASI experiment [72], the CAPMAP results [73], the (almost) three year results of CBI 
[HI and the preliminary results of QUAD [73]. In the present version of the code the polarization 
correlations and cross-correlations (i.e EE and TE power spectra) can be explicitly computed since 
consistent initial conditions have been given for the whole Boltzmann hierarchy. The magnetic field 
can have two distinct effects on the CMB polarization. Since gravitating magnetic fields modify the 
structure of the adiabatic mode and of the evolution of the baryon-photon fluid, the TE and EE 
angular power spectra will be different. The second effect would be due to the presence of a Faraday 
rotation term which would couple the evolution equations for the two brightness perturbations which 
are sensitive to polarization. In the language of Eqs. (12.811) and (12.821) this term would couple the U 
and Q Stokes parameters producing, ultimately, a rotation of the polarization plane of the CMB. The 
Faraday coupling can be easily included if the magnetic field is uniform [75]. In the case of stochastic 
magnetic field this calculation has never been done. There are certainly semi-analytical attempts in 
this direction (see, for instance, [77t [78l 179] ). However, the main problem with these calculations is 
that they assume that the magnetic fields only rotate the polarization without entering in any other 
place of the evolution equations. There are, on the contrary, reasons to believe that, for a stochastic 
field, the two effects can be equally important^. In this discussion we will not include the Faraday 
rotation term by a uniform field since this would break explicitly the spatial isotropy and concentrate 
on the numerical solution when the magnetic fields are consistently introduced in the initial conditions 
and in all the other evolution equations, as it was done for the TT correlations in the previous section. 

In Fig. [HI the TE correlations are illustrated for the case of blue magnetic spectral indices 1 < n-Q < 
5/2. The parameters used in the calculation are exactly the ones employed in Fig. [H If the magnetic 
field is of the order of 0.1 nG the magnetized TE correlations cannot be distinguished from the three 

^^In nearly all Faraday rotation studies it is assumed that the description of the plasma is given in terms of a single 
fluid. This assumption is not correct |80| . Indeed Faraday rotation requires necessarily a kinetic treatment (or at least 
a two fluid treatment). We are investigating the possibility within our numerical approach [50j . 
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Figure 8: The height of the first and second peaks for different values of the magnetic field intensity 
as a function of the magnetic spectral index. The other cosmological parameters have been fixed as in 
Eq. dHD. 



year best fit of the WMAP data. Unlike the case of the temperature autocorrelations (where the 
position of the Doppler peak cannot be moved by a stochastic magnetic field) there is an observable 
shift of the second and third (correlation) peaks of the TE spectra. This distortion also entails a shift 
of the position of the corresponding peaks. A similar effect can be observed in the magnetized EE 
correlations which are reported in Fig. [71 Also in this case the peaks are raised and partially shifted. 
Figures Eland [7] show, a posteriori, that the magnetic fields also affect the polarization observables even 
without a Faraday rotation term. This observation supports our previous statements. The physical 
reason of the obtained result can be understood very simply. To zeroth-order in the tight-coupling 
expansion, the magnetic field affects the dipole of the brightness perturbation for the intensity. Always 
to zeroth order, this contribution is reflected in a further source term for the monopole. But both the 
TE and EE power spectra arise to first-order in the tight-coupling expansion and are proportional to 
the first-order dipole through a term which is, up to a numerical factor, k/e' [32]. This shows why we 
also get an effect on the polarization observables even if the Faraday rotation term is absent. 

The results obtained so far show that it is possible to obtain accurate estimates of the temperature 
autocorrelations and of the polarization correlations also in the presence of a magnetized background. 
Conventional CMB calculations have a high level of accuracy and this is due, both, to the precise 
understanding of the initial conditions and to the thorough comprehension of the dynamics. At the 
moment, it is possible to achieve the same level of accuracy also when a magnetized background is 
included. 



6 Waiting for Planck 

The Planck explorer satellite [22] will provide high precision measurements of the cosmic microwave 
background. In view of this exciting time the quantitative trends illustrated in Sections jH and 
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Figure 9: The ratio of the heights of the second to the first peak (plot at the left) and of the heights 
of the third to the first peak (plot at the right). The parameters are fixed as in Fig. [Hi 

will now be scrutinized in more depth. The Doppler peaks are both distorted and increased. For 
2 < £ < 2500 the extrapolated best fit to the WMAP data alone predicts, at different locations, 7 
acoustic peaks. The heights and shapes of the 7 peaks have been monitored for different values of the 
magnetic field intensity and of the magnetic spectral index. For reasons of space we will just focus on 
the first and second and on the sixth and seventh. This will suffice for the quantitative trend we wish 
to illustrate. 

In Fig. [8] the height of the first peak is reported for different values of the regularized magnetic 
field, as a function of the magnetic spectral index. Note that in all the figures from Fig. [8] to [12] 
the title of each plot labels the peak. So, for instance, the notation "1. Doppler peak" means that the 
corresponding plot refers to the first Doppler peak. 

The maximal increase of the acoustic peaks always arises for intermediate spectral tilts. The TT 
correlations are not only shifted upwards but they are also distorted: this is evident from Fig. M 
where we illustrate, respectively, the ratio of the second peak to the first (plot at the left) and the 
ratio of the third peak to the first. From the left plot of Fig. [9] it is clear that an increase in the 
spectral index entails, for ub < 1-8, a decrease of the height of the second peak in comparison with the 
first. When ub > 1-8 the opposite behaviour is observed. For a given value of the magnetic field the 
correlated distortion of the first peak is an effect of the order of few percent (as it can be argued from 
the corresponding plots). As we reach into the region i > 1500 the effect becomes more pronounced 
especially for the seventh peak. While the first peak is not sensitive to a nG magnetic field, the seventh 
peak can be a reasonable indicator of the presence of large-scale magnetic fields in the nG range. This 
aspect is illustrated in Fig. [TDjin the cases of the sixth and seventh peaks. From the plot at the right, 
for instance, a 15 nG field has a 10 percent effect on the shift. The precise value depends on the 
spectral index, as it can be argued from Fig. [TOl 

The observed correlated distortion is also visible in the case of the TE angular power spectra. In 
Fig. [11] the first anticorrelation peak (plot at the left) is compared with the first correlation peak (plot 
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Figure 10: The heights of the sixth (plot at the left) and of the seventh (plot at the right) Doppler 
peaks is illustrated for different values of the regularized magnetic field as a function of the spectral 
index. 

at the right). The nominal value of the first correlation peak (appearing in the title of the plot) is the 
one stemming from the three year best fit to the WMAP data alone. 

In Section H] we pointed out that a slight increase in the CDM fraction could be compensated by 
the presence of a magnetic field. This kind of potential degeneracies can only be throughly discussed 
in the framework of a general parameter estimation which also includes, to begin with, the magnetic 
field parameters. This analysis is beyond the scope of this paper, however, it is useful to investigate in 
an eclectic perspective, also other potentially interesting degeneracies which can be only assessed (or 
even partially resolved) in more systematic approaches. 

One of these potential degeneracies involves the optical depth to reionization, i.e. r. The increase 
of r in a model without magnetic field yields a lower height of the Doppler peak. The values of the 
best fit model for WMAP data alone are used (see Eq. (14. ip ) apart from the optical depth r which is 
assumed to be varying between 0.09 and 0.105. In Figure [12] the height of the first acoustic peak is 
shown for different values of the magnetic spectral index ns for different values of the magnetic field 
strength and of the optical depth r. 

The dependence on the magnetic spectral index becomes more important as -Bl increases. While 
for i?L = 2 nG the height of the acoustic peak is basically independent of the spectral index, there 
is, for i?L = 16 nG, a nontrivial functional relation between the height of the peak and the magnetic 
spectral index. The variations in r can be partially compensated by the addition of the magnetic field. 
Note, indeed, that the full line denote the experimental indetermination in the position of the peak. 

In the near future the Planck explorer mission with its low frequency and high frequency instru- 
ments will be able to probe with higher accuracy the temperature autocorrelations, the polarization 
autocorrelations and cross- correlations. At the moment the putative Planck sensitivity can only be 
inferred from the last version of the Planck blue book p2]. The nominal sensitivity might not be the 
one effectively achieved by the instruments. Given the present specifications of the instruments, is not 
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Figure 11: The height of the first anticorrelation and correlation peaks is illustrated, respectively, in 
the left and in the right plots. 



clear in our case what are the best frequency channels to be combined in order to be mostly sensitive 
to the features introduced in the angular power spectra by the large-scale magnetic fields. In spite of 
this we find it interesting to elaborate on the possible implications of our endeavors for the Planck 
measurements at high multipoles (i.e. i > 1500). 

To compare the nominal Planck sensitivity with the situation where magnetic fields are consistently 
included in a ACDM paradigm (with no tensors) we find it useful to adopt the following measure: 

^(TT) ^ \Ci''^\Bun^)-Cr\B^ = 0)\ 

""'^ Cf^^(SL = 0) 

In Eq. (16. ip C^^^\Bi^ = 0) is computed from the three year best fit to the WMAP data alone; 

(Slj'^b) is the TT correlation but computed with a magnetic field of regularized intensity i?L 
and characterized by a spectral index n^. In different frameworks a similar estimator has been also 
employed [231 [2l]. 

In Fig. [13] the quantity defined in Eq. (16. ip is illustrated for different values of the magnetic field 
intensity and of the spectral index. Note that q estimates the difference induced by the presence 
of the magnetic field on the extrapolated three year WMAP best fit which can be used to deduce the 
nominal sensitivity of Planck for different regions in the multipole space. Adopting the three year 
WMAP best fit as fiducial model, the l-cx errors can be inferred following the standard analysis also 
thoroughly reviewed in the Planck blue book (see [22]) 

WtW! = j;4"e-'<'+'>'';~.'. (6.3) 

C 

In Eqs. (I6.2p - fl6.3p various assumptions should be made as far as the sky coverage and the relevant 
frequency channels are concerned. For this reason what we are presenting here are just preliminary 
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Figure 12: The height of the Doppler peak is presented as a function of the magnetic spectral index ub 
for different values of the magnetic field strength and the optical depth. The horizontal lines indicate 
the observational bounds on the Doppler peak from WMAP3. 

indications of what could be the trend of the Planck accuracy on the basis of the figures customarily 
employed by the Planck team [22] . Needless to say that the present estimate can be made more realistic 
once the effective Planck sensitivity will be available. Thus, it will be assumed that /sky — 0.65 
corresponding to ±20° galactic cut. Furthermore the three lowest frequency channels of the high 
frequency instrument (i.e. 100 GHz, 143 GHz and 217 GHz) are combined and the sum appearing in 
Eq. (16. 3p then extends over these three channels. To make explicit the sum it should be noted that 
= ('^pt'^^fwhm)"^ is the sensitivity per resolution element 'J^fwhm ^ ''^fwhm- The quantities a^rj, and 
''^FWHM change for each of the three aforementioned channels. For a^j, and "(^fwhm ^^e values reported 
in the Planck blue book [22] have been selected. With these specifications in mind, the last quantity 
to be defined is ^bea,m ~ V8 In 2 / ^fwhm which measures the resolution of the Gaussian beam. 

In Fig. [13] with the full curve the 1-a error is reported. If the estimator leads to a value which is 
larger that the foreseen sensitivity it will be possible to make observational distinction between the 
magnetized model and the extrapolated three year best fit. In spite of the intrinsic uncertainty on the 
actual sensitivities of the instrument Fig. [13] is eloquent enough and then offers encouraging prospects 
for the region of high multipoles. 

7 Concluding remarks 

Large-scale magnetic fields are an observed component of the present Universe. Why are they present 
at the fiG level in galaxies with different morphologies and different evolutionary histories? Why are 
they present inside rich clusters? When did they originate? Are they a cosmic relic in the same way 
as the adiabatic mode of curvature perturbations has a primeval origin? 

It is difficult to answer these rather motivated questions just by building new models or by empha- 
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Figure 13: The estimator of Eq. (16. ip is illustrated in combination with the foreseen Planck sensitivity 
(full curve). 

sizing single (potentially interesting) effects. It is even more difficult to run complicated simulations 
trying to reproduce large-scale magnetic fields without knowing the initial conditions to be imposed at 
early times, i.e. after decoupling. What we need is a systematic scrutiny of pre-decouphng physics to 
answer more modest but necessary questions. The most urgent one concerns the effect of large-scale 
magnetic fields on CMB observables. The combined numerical and analytical tools presented in this 
paper are a promising step along this direction. Needless to say that we plan to scrutinize more deeply 
all the manifold and exciting implications of our approach. 
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A From synchronous to longitudinal gauges 



The treatment employed in the present analysis is the one stemming from the synchronous gauge 
description. We cross-checked all our results within the longitudinal approach. This cross-check is, 
under certain circumstances, mandatory. Indeed, in the synchronous gauge the freedom of selecting 
the coordinate system is not completely fixed. It is therefore important to have at hand a gauge 
description where the gauge parameters are completely fixed. 

The strategy we followed has been often to derive the same quantity in the two different gauges 
and then compare the results by transforming the obtained expressions from one gauge to the other. 
For the effectiveness of this approach the longitudinal gauge is not essential. Rather we should say 
that the only essential requirement is a gauge where the freedom of selecting the coordinate system 
is completely fixed. So, for instance, the uniform curvature gauge would work equally well for this 
purpose [m |82] . 

In the longitudinal gauge the metric of Eq. (12. ip is perturbed in such a way that non-vanishing 
entries of the first-order metric are 

5^900 = 2a^(r)0(/i;, r), S^gij = 2a^(r)V^(A;, T)5ij. (A.l) 

The difference between the longitudinal and the synchronous coordinate systems is evident by compar- 
ing Eq. (lA.ip with Eq. (12. 7p . Following standard techniques we can find the precise relation between 
the longitudinal and the synchronous degrees of freedom: 

Hk, r) = -^{[h{k, r) + Qi{k, r)]" + UHk, r) + Qi{k, r)]'}, 
^(fc, r) = -i{k, r) + ^^[h{k, r) + Qi{k, r)]', 
-5{k, r) = 6{k, t) + ^'^^^jj~^\ h{k, t) + Qiik, r)]' 

e{k, t) = eik, t) - ^[h{k, r) + 6e(fc, r)]'. (A.2) 

The barred quantities (i.e. 6 and 6) are defined in the longitudinal gauge and w is the barotropic index 
of the corresponding species. Similarly the transformation for 6 holds for a generic peculiar velocity. 
The inverse transformations can be also obtained and they are: 

e(fc, r) = -^(A;, r) - - r a(r')</)(A:, r')rfr', 

a J 

h{k, t) = 6^{k, r) + 6- / a{T')(P{k, T')dT' - 2^ / — — / a(r")</>(fc, r")rfr", 

a J J fl(T ) J 

6{k, t) = 5(fc, r) + ^^^"^ ^ r a{T')<p{k, T')dT' 

a J 

1 2 

eik, t) = eik, r) r air'Uik, T')dT' . (A.3) 

a J 

The integrals appearing in Eq. (lA.Sp for the expressions of e and h imply two integration constants 
which can be space dependent and which are fixed by demanding that = and that h has no 
constant mode. 
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The solution for the magnetized adiabatic mode will be, in the longitudinal gauge, 



20C{k) _ ^RjlAa^ik) - RMk)] 



i'*{k) 



AR^ + 15 
/8/?^ + 20 
\AR, + 15 



C{k) + 



AR^ + 15 
R^[AaB{k,T) - R,^^{k)] 



AR^ + 15 



R^ 



iP,{k)=[l + -R, U,{k) + ^[AaB{k) - R,n^{k)] 



d^{k,T) 


= -2(f),{k) - R^n^ik), 




5u{k) = 


~2Mk) - R^^Bik), 




Scik) = 


3 3 

- -Rj^lB{k), 




S^{k) = 


3 3 




(Ty{k,T) 


Ry 

Ry 




0^h{k,T 


k^T 

2 


\Mk) + 2Mk)\ 






2 


Mk) "^^^2"^^^+'^'"^^ 


0c{k,T) 


k^r 
= 2 







(A.4) 

It can be easily checked that the solution (1A.5P is a solution of the full system written in the longitudinal 
frame. 

There is an important point to be borne in mind when setting initial conditions. We set initial 
conditions deep in the radiation epoch. Now, the constant C = C{k) that appears in the synchronous 
description can be actually related with curvature perturbations on comoving orthogonal hypersurfaces 
TZ. In the longitudinal gauge we do know that, deep in the radiation epoch. 



n{k) = -ip 



-i'*{k) 



<p*{k) 



-2C{k) 



(A.5) 



In the longitudinal gauge we can also easily express the variable C, introduced in Eq. f l2.73p . In terms 
of the longitudinal degrees of freedom 



c 



5^Pt + 5pB 



(A.6) 



where 5^pt is the total density fluctuation in the longitudinal gauge. In the longitudinal gauge the 
Hamiltonian constraint reads 

- miUct) + tl)') = AirGa^iS^pt + Sps). (A.7) 
Using Eqs. ( ]A.5p and ( ]A.6p into Eq. ( 1A.7I) the Hamiltonian constraint can be written as 



C = n 



127rGa2(pt + Pt) 



(A.8) 



where Eqs. (12.31) and (12.41) have been used. Equation ( lA.SI) has been quoted and independently 
obtained in Eq. (12.751) 
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